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Balloon-Born^  High-Altitude  Gravimetry 
The  Flight  6f  DUCKY  la 
fll  October  1983) 


1.  INTRODUCTION 

Gravity  field  values  at  high  altitudes,  between  altitudes  fjr  aerial  surveys 
and  satellite  orbits,  (e.g. ,  30  km),  are  normally  estimated  from  upward  continua¬ 
tion  of  surface  measurements,  downward  continuation  of  satellite  measurements  as 
computed  from  orbital  perturbations,  and/or  geoid  models.  These  techniques  are 
well-developed  with  generally  accepted  results,  but  are  subject  to  some  limitations. 
Upward  continuation  depends  upon  the  quality  and  distribution  of  surface  data, 
usually  nonuniformly  spaced,  and  taken  from  different  surveys;  thus,  gaps  and 
uneven  spatial  distribution  may  result  in  inaccurately  upward  continued  estimates. 
With  increasing  altitude,  short  wavelength  information  on  crustal  structure  is 
attenuated,  not  recoverable  from  downward  continuation  from  satellite  altitudes. 
Geoid  models  tend  to  emphasize  wavelengths  longer  than  about  30  km,  and  are 
well-known  only  over  oceanic  areas. 

Verifying  gravity  data  bases  and  models  requires  establishing  the  validity  of 
the  models  at  locations  where  measurements  have  not  been  made.  While  such 
models  currently  can  be  tested  effectively  at  ground  level,  their  validation  at 
altitude  awaits  the  development  of  a  suitable  approach.  The  Air  Force  Geophysics 
Laboratory  (AFGL)  is  developing  a  program  to  verify  gravity  model  estimates  by 
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measuring  gravity  directly  using  high-altitude  balloons.  The  basic  concept  is 
shown  in  Figure  1.1.  A  gravimeter  package  suspended  beneath  a  balloon  is  in  a 
dynamic  and  largely  unpredictable  environment,  sensing  not  only  the  gravitational 
acceleration  of  the  earth  but  also  all  accelerations  due  to  the  motions  of  the  balloon 
system.  An  exaggerated  sketch  of  kinematic  balloon  motion  is  shown  in  Figure  1.2. 
For  a  specified  time  interval  (e.g. ,  1  second)  during  which  a  measurement  is 
made,  the  variation  in  balloon  accelerations  is  expected  to  be  significantly  greater 
than  the  variation  expected  in  the  Earth's  field.  Therefore,  additional  instrumen¬ 
tation  is  required  to  measure  as  many  balloon  motions  as  possible,  such  as  rotation, 
bobbing,  and  swaying.  As  all  such  ancillary  sensors  are  dependent  on  the  local 
inertial  frame,  gravitational  acceleration  cannot  be  separated  from  vertical  balloon 
accelerations  without  additional  data  acquired  independently  of  the  balloon's  ref¬ 
erence  frame.  These  independent  data  are  extracted  from  balloon  tracking  (Figure 
1.3)  which  must  accomplish  three  objectives:  1)  the  measurement  of  the  gravi¬ 
meter  package  accelerations  (especially  the  vertical)  referred  to  a  ground-based 
coordinate  system;  2)  the  measurement  of  velocity  for  estimation  of  the  Eotvos 
effect;  and  3)  the  measurement  of  the  gravimeter  position,  which  is  used  as  an 
input  to  the  gravity  model.  Combining  balloon  data  with  tracking  data  allows  the 

separation  of  balloon-induced  accelerations  from  gravitational  accelerations. 

~3  2 

The  long-term  goal  is  to  determine  gravity  to  1  mCal  (10  cm /sec  ),  and  in 
support,  obtain  instrument  motion  data  (worst-case  limits)  of  accelerations  to 
1  mGal,  velocity  to  5  cm /sec,  and  position  to  3  m  in  all  three  orthogonal 
coordinates. 

The  balloon  launch  was  siet  to  coincide  with  the  lowest  seasonal  wind  velocities 
(Figure  1.4)  over  Holloman  Air  Force  Base  where  AFGL  has  its  permanent  balloon- 
launch  facility.  The  mild  wind  velocities  are  desired  to  provide  the  most  benign 
environment  possible  during  the  testing  phase  so  that  the  analysis  effort  may  focus 
on  the  performance  of  the  gravimeter  system.  Balloons  have  been  flown  up  to 
150,000  ft  where  the  eir  pressure  is  less  than  1  mm  of  Hg.  Temperatures  can 
range  from  -70°  C  at  night  to  almost  40°  C  on  exposed  parts  of  the  gondola's 
aluminum  framework  in  the  daytime.  The  instrument  packages  had  to  withstand 
the  shock  and  vibration  associated  with  both  launch  and  landing.  Even  with  a 
parachute  and  a  crush  pad  on  the  gondola,  accelerations  up  to  8  g  may  be 
expected  on  landing.  Ideally,  the  instrument  packages  would  be  self-contained 
except  for  telemetry  and  batteries. 

The  instruments  were  also  operated  in  a  strepdown  mode;  i.e.,  the  instrument 
axes  are  not  fixed  with  respect  to  inertial  space  liu  are  attached  to  the  instrument 
package  itself.  Since  the  packages  are  rigidly  fastened  to  the  balloon  gondola, 
the  instrument  axes  correspond  to  a  local  coordinate  system  in  the  gondola. 
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GRAVITATIONAL 

ACCELERATION 


Figure  1.2  A  Sketch  Showing  Exaggerated  Motions  of  the  Instrument  Package. 
Each  of  these  motions  contributes  to  a  vertical  acceleration  term;  thus,  con¬ 
fusing  the  gravimeter  data 
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MEAN  EAST- WEST  WIND  COMPONENTS  (knots) 
Holloman  AFB 


2.  EXPERIMENTAL  DESIGN 


2.1  Gondola  Design  Concept 

General  requirements  for  a  gravity  measurement  flight  system  began  to  devel¬ 
op  in  March  1983.  A  target  flight  data  for  the  completed  system  was  set  for 
October  1983.  This  relatively  short  lead  time  appeared  reasonable  because  devel¬ 
opment  of  a  balloon-borne  motions-sensor  package  had  already  been  started  under 
an  AFGL  contract  with  the  University  of  Toronto.  This  package  would  contain 
accelerometers  and  strapdown  rate  gyros  aligned  in  three  orthogonal  axes  for 
determining  altitude  and  motions  of  the  payload. 

The  scientific  goal  was  to  resolve  accelerations  in  the  X,  Y,  and  Z  planes  to 
one  part  in  10  .  This  meant  every  effort  had  to  be  made  to  minimize  extraneous 
motion  and  to  attempt  to  achieve  a  payload  configuration  that  would  hang  beneath 
the  balloon  as  close  to  vertical  as  possible.  In  essence,  we  should  try  to  produce 
a  true  plumb  bob  when  airborne.  We  also  selected  a  symmetrical  shape  to  reduce 
the  effect  of  wind  and  turbulence. 

Attachments  were  designed  to  limit  the  tendency  of  a  lengthy  flight  train 
(balloon-recovery  parachute  and  gondola)  to  wind  up  and  unwind  as  the  system 
ascends  to  float  altitude.  During  ascent,  the  balloon  changes  shape  as  it  rises 
and  imparts  rotation  to  the  system,  storing  this  rotational  energy  in  the  in-line, 
elastic  nylon  recovery  parachute.  The  mechanical  aspects  of  the  flight-system 
design  thus  became  paramount  considerations  to  minimize  undue  linear  motions  and 
rotation.  Essentially,  we  limited  the  degrees  of  freedom  in  the  rigging,  and  de¬ 
signed  for  maintenance  of  balance  and  verticality  during  operation  of  any  on-board 
control  devices  such  as  ballast  dispensing  and  valving  off  helium  lifting  gas. 

Soon  other  mechanical  problems  surfaced,  due  to  the  need  to  physically  and 
electrically  integrate  the  tracking  beacons,  the  radar  transponder,  the  recovery 
beacon,  the  two  down-link  data  transmitters,  the  six  antennas,  and  a  full-time 
command-control  capability,  and  to  provide  internal  and  external  power  system 
compatibility.  Temperature  controls  for  the  prime  sensing  devices  were  required 
to  minimize  drift  and  subsequent  error  in  measurement. 

The  melding  of  several  independent,  though  closely  interrelated,  subsystems 
necessitated  devising  a  gondola  system  that  could  provide  sectional  disassembly  and 
operational  independence  for  service  and  individual  testing/calibration,  yet  permit 
flexibility  for  modification.  The  size  and  form  factor  of  the  gondola  was  also 
limited  by  the  launching  method  to  be  employed  for  the  balloon  system.  A  dynamic 
launch  method,  with  the  gondola  suspended  from  the  boom  of  a  mobile  P&H  30-ton 


crane,  limited  the  vertical  size  that  could  be  handled  safety.  Using  a  balloon  of 

2.1  million  cubic-feet  volume,  the  payload  weight  had  to  remain  under  2,500  lbs 
including  445  lbs  of  pourable  glass-bead  ballast,  in  order  to  reach  the  desired 
float  altitude  and  achieve  the  desired  flight  profile  of  two  sampling  altitudes. 

The  most  suitable  shape  for.  meeting  ail  requirements  was  a  long  cylinder  that 
consisted  of  five  sections  bolted  together  oh  internal  flanges.  All  sections  were 
34  in,  in  diameter.  Section  heights  varied  from  23.5  In.  to  33  in.  for  a  total 
stacked  height  of  il  ft  3  in.  A  four-point  suspension,  using  all-stainless  cable 
(wire  rope)  with  adjustable,  fine  thread  turnbuckies  in  each  leg,  terminated  into  a 
singlepoint  fixture  for  attachment  to  the  launch  crane's  tri-plate  release  device. 

The  parachute  was  specifically  oversized  to  siow  the  impact  velocity  to  15  to 
20  ft/sec,  and  a  crush  pad  of  stepped,  corrugated  cardboard  16  in.  thick  was 
attached  at  the  base  of  the  cylindrical  gondola.  A  roll  bar  timpact-absorbing 
brush  ring)  7  ft  in  diameter,  concentric  with  the  cylinder  axis,  assured  that  with 
anticipated  horizontal  velocities  of  30  to  50  ft/sec,  the  gondola  would  lie  over 
safeiy  with  minimal  physicdl  damage  upon  ground  impact.  The  base  of  the  crush 
pad  carried  impact  switches  designed  to  release  one  set  of  parachute  suspehsipn- 
harnesS  risers  to  collapse  the  chute  and  prevent  the  possibility  of  dragging  the 
payload  across  the  ground  after  impact. 

the  exterior  of  the  cylindrical  gondola  was  skinned  with  aluminum.  When 
flight-ready,  it  was  wrapped  with  2-in.  thick  fiberglass  "Certalnteed"  pipe  in¬ 
sulation,  with  an  All  Service  Jacket  (ASJ)  vapor  barrier  of  aluminum  foil,  and 
glass-reinforced,  whitecoated  Kraft  laminate,  which  served  as  a  reflective,  white 
outer  surface,  finished  with  superior  compression  resistance  for  added  protection. 

2.2  Flight-Control  and  Data  Systems 

The  design  of  the  flight-control  system  followed  standard  AFGL  practice  of 
providing  primary  and  backup  systems  to  ensure  positive  experiment  control  and 
redundant  flight-safety  control  via  separate  power  systems,  in  this,  case,  because 
of  the  stacked  sectional  design  of  the  gondola,  the  independent  primary  and  back¬ 
up  systems  were  located  in  the  same  gondola  section.  One  section  was  to  contain 
ail  necessary  battery  power.  This  was  achieved  by  using  two  "freeze-  pack" 
thermal  containers,  each  housing  fourteen  80  Ah  BB622/U  batteries,  configured  in 
such  a  fashion  as  io  provide  isolated  power  sources  for  the  various  ilisiruments 
and  sensors  aboard  the  gondola. 

There  were  five  separate  internal  power  sources,  plus  one  with  the  capability 
of  continuous  external  power.  The  power  system  weighed  approximately  322  lbs. 
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Each  power  source  could  be  powered  up  or  down  by  radio  command  of  either  .the 
primary  or  backup  digital  command  systems. 

A  single. instrument  frame,  17  x  17  x  10}  in.,  housed  all  flight-control  elec¬ 
tronics,  two  digital  command  units,  housekeeping  sensors  and  monitors,  power 
control  interface  units,  and  one  PCM  (pulse-code  modulation)  data  encoder.  This 
instrument  frame  was  housed  in  a  laboratory-grade,  low-temperature-insulated 
container  (bio  specimen  "freeze-pack”)  for  environmental  protection  and  tempera¬ 
ture  control.  This  method  of  packaging  permits  very  high  flexibility  as  the  final 
system  configuration  is  composed  of  an  assembly  of  modular  items,  all  intercon¬ 
nected  by  a  master  harness  terminated  in  a  master  patcn-panei-type  connector 
interface  for  inputs  and  outputs.  ..  - 

The  gondola  had  two  telemetry  down  links  and  two  radio-command  up  links. 
Both  telemetry  links  were  continuous,  and  employed  standard  IRIG  (Inter  Range 
Instrumentation  Group)  PCM/FM  data  format.  The  housekeeping  and  motion- 
sensing  data  were  sampled  at  128  kilobits  per  second,  using  ten  bits  per  word, 
most  significant  bit  sent  first;  there  were  47  analog  words  (channels),  and  three 
digital  words  with  two  synchronizing  words  without  parity  check,  included  per 
data  frame.  The  output  code  was  transmitted  Bi-Phase  FM  on  a  500  kHz  modulated, 

2.2  GHz  carrier  with  2  w  RF  output. 

The  science  data  from  the  Vibrating  String  Accelerometer  (the  prime  gravity 
sensor),  with  its  own  housekeeping  and  monitoring  information,  was  transmitted  on 
a  separate  2.2  GHz  PCM  channel  due  to  the  nature  of  the  sensor  output  and 
sampling  rate  desired.  For  this  we  utilized  a  19.2  kilobaud  bi-phase  encoder  of  30 
words  per  frame  (27  analog  and  3  digital).  This  was  also  a  standard  IRIG  format, 
10  bits  per  word. 

On  the  ground,  both  data  streams  were  recorded  on  wideband  FM  analog  tape, 
and  simultaneously  stripped  out  for  display  in  real  time,  on  digital  meters  or 
graphic  chart  recorders.  Thirty-six  strip  chart  data  channels  were  displayed  in 
real  time  during  the  flight  to  assess  program  progress.  A  complete  listing  of  the 
data  channel  assignments  appears  in  Appendix  B. 

2.3  Scnsorg 

2.3.1  VIBRATING  STRING  ACCELEROMETER  (VSA) 

The  Vibrating  String  Accelerometer  (VSA)1  is  a  single-axis  accelerometer 
originally  designed  for  the  Atlas  missile  navigation  system  in  the  1960s.  When  the 
VSA  was'  deemed  obsolete  for  the  Atlas,  the  Woods  Hole  Oceanographic  Institution 


1.  Bowin,  Carl  O.,  Scheer,  Eddie,  and  Goldsborough,  Rob  (1984)  Technical 
Report  -  Balloon  Environmental  Test  Report,  AFGL  Contractor  Report. 
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(WHOI),  Woods  lip|e»  MA  acquired  a  collection  of  YSAs  from  government  surplus 

stock.  WHO?  developed  and  byilt  a  shlpborne  gravimeter  system  based  op  the 

'234” . .  . .  ”  "■  . *"  '  ■  ’  •  *  . . 

VS  A 1  ’  and  successfully  used  it  in  geophysical  research  over  the  subsequent 

years.  When  the  balloon-gravity  program  pame  Into  being,  WHOI  proposed  ?o  build 

a  graylmeter/inertial  system  fpr  AFPL,  In  fylarph  1983 ,  WHOI  came  under  contract 

to  AFGL  and  successfully  built  0  gravimeter  system  ip  ?|me  for  the  Planned 

September  1883  launch, 

2.9.1.1  Mephapical  hayoyt  of  tit?  System 

Six  basic  modules  were  developed  at  WHQI  apd  placed  aboard  til?  gondola: 

1)  The  VS  A  and  Oven  System:  a  fcan'1  (VSA  two-staqe  ovpn)  hpusing  th?  VS  A  and 
its  temperature-control  components  and  circuitry;  2)  The  Processing  System:  a 
cardcage  holding  six  circuit  cards;  3)  The  Oscillator:  a  shielded  ’tpyho?  holding 
two  mpre  cards:  a  10  MHz  crystal  oscillator  apd  associated  oscillator  control  cir¬ 
cuitry;  4)  The  Encoder:  (obtained  from  AFGL)  that  relayed  acquired  Y^A  data  to 
ground;  5)  The  Heat  Sink:  e  large  heat  ?lnH  upon  which  wero  mounted  both  a 
power  supply  that  converted  Rower  from  the  AFfih  batteries  for  use  in  ysa  cir¬ 
cuitry,  end  a  set  Of  power  transistors  used  to  drive  YSA  heater  blankets;  and  6) 

The  Front  Panel:  ysed  for  monitoring  and  adjusting  vsa  component®  while  the 
system  was  on  the  ground! 

in  addition,  a  Micromatlon  mleroeompHtert  a  Kennedy  tap?  drive*  an  eight-digit 
display,  and  assorted  power  syppiies  and  test  equipment  were  supplied  by  WHPt 
for  testing,  and  fnr  acquisition  pf  data  taken  on  th?  ground  prior  to  the  experiment* 

2.34.2  vsa  System  Pespription 

Figure  24  is  a  block  diagram  schematically  illustrating  the  basic  components 
and  interconnections  amongst  the  VSA  system  components. 

VSA  and  Oven  System 

The  principal  components  of  the  cap  are:  1)  VSA  Sensor;  The  VSA  is  a 
device  consisting  of  a  mass  suspended  by  two  m?tal  band?  (strings)  under  tension 

(figure  2.2).  The  tension  on,  and  the  length  of,  eaoh  string  Is  set  so  th?  r?®o~ 
nant  frequency  is  about  j  kfiz  at  a  very  high  9  (vald?  not  available)*  PH?  to  the 
high  9,  the  string®  will  oscillate  (at  rasonsnee)  vary  easily  t  and  amhlant  noise 


2.  Bowin,  Pari,  wing,  Phsrles  fi.»  and  Aidriph,  Thomas  p;  (19§0)  Test  of  the 

MIT  vibrating  strjng  gravimeter i  J.  Geophys.  Res.  74,  12:3278-3280. 

3.  Bowin,  par?  q.,  Fpiinsbee,  Allln>  and  Aidrioh,  Thomas  P.  (1970)  Test  of 

digital  YSA  spa  gravity  meter  and  comparison  with  hsposte  apd  Bpmberg 
gyrpstabilized  gravity  meter  (Abstract),’ Trans.  Am.  fiedphys.  Union  5f, 
4:2gi,  ..........  •  ‘  - . - 

4.  Bowin,  p.,  Aidriph,  t*  c,,  and  Foiinabee',  $.  (197?)  ysa  gravity  meter 

eyet'em:  tests  and  recent  developments,  J.  Geophys.  Res.  77: 2018-2033 . 


VSA  SENSOR 


f, ,  f2  ~  4  KHz  f,  -  f2  ~  60  Hz  2 

AT  9.8  m/*ec 


Figure  2.2  A  Sketch  of  the  VSA  Sensor.  Two  strings  under  tension  are  separated 
by  a  mass.  A  force  acting  on  the  mass  along  with  the  strings  will  cause  different 
tensions  on  each  string;  thus,  each  string  will  have  a  different  resonant  frequency. 
The  VSA  system  measures  the  acceleration  from  the  beat  frequency  of  the  two 
strings.  The  resonant  frequency  of  each  string  is  about  4  xHz;  the  beat  frequency 
at  1  g  is  about  60  Hz 
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is 'enough  to  drive  the  oscillations .  The  surrounding  magnets -sustain- the  oscil¬ 
lations.  The  oscillation  creates  an  electromagnetic  force  (EME)  due- to.  the  ambient 
magnetic  field.  The.  EMF  is  sensed,  amplified,  and  fed  back  to  the  strings.  When 
an  acceleration  is  sensed  in  the  sensitive  axis  (along  the  strings),,  the  force  acting 
on  the  mass  results  in  different  tensions  acting  on  the  two  strings;  The  string  in 
the  direction  of  the  acceleration  is  at  a  lower  tension  than  the  opposite  string.  At 
1  g,  the  difference  in  resonant  frequencies  is  about  60  to  64  Hs.  Thus,  when  the 
two  frequencies  are  mixed,  three  narrow  bands  of  frequencies  exist:  65  Hs,  4  kHs, 
and  8  kHs.  When  low-pass  filtered,  the  final  VS  A  output  frequency  is  about  65  Hs, 
and  it  is  this  frequency  that  is  a  measure  of  the  along-axls  acceleration.  For  this 
flight,  we  chose  to  sample  the  VS  A  output  every  10  cycles,  giving  an  effective 
sampling  rate  of  about.  6,5  Hz.  The  circuitry  was  designed  to  yield  either  a 
frequency  count  or  the  period,  the  latter  being  the  choice  for  the  flight; 

2)  Osc/Amp:  An  Oscillator/Amplifier  circuit  that  effectively  sustains  oscillations  in 
one  of  the  two  strings  of  the  sensor  by  positive  feedback.  The  output  of  the 
circuit  is  a  sinusoidal  signal  with  a  nominal  frequency  of  4.5  kHs.  There  are  two 
of  these  circuits  in  the  can,  and  two  strings  in  the  sensor.  The  difference  in  the 
frequency  of  the  two  resulting  signals  is  directly  proportional  to  the  acceleration 
and/or  the  component  of  gravity  experienced  along  the  input  axis  of  the  sensor: 
a=K*fd  +  B,  where  a  is  acceleration  in  milligals,  K  is  the  scale  factor,  fd  is  the 
difference  frequency  in  Hz,  and  B  is  a  bias  term  in  milligals.  Prior  to  the  experi¬ 
ment,  a  scale  factor  determination  was  performed,  and  the  following  values 
obtained: 

K  =  15218.29659  (milligals/Hz)  (2.1) 

B  =  4980.720047  (milligals)  (2.2) 


A  rough  value  for  fd  at  1  Q  acceleration  *8  64.1  Hz.  3)  Inner  and  Outer  Ovens: 

A  two-stage  oven  is  used  to  maintain  the  temperature  of  the  sensor;  the  "outer" 
housing  the  "inner."  Each  oven  contains  two  thermistors  that  sense  temperature 
within  the  controlled  volume:  a  "monitor"  thermistor  that  is  used  solely  for  tem¬ 
perature  observation,  and  a  "control"  thermistor  that  is  connected  in  a  bridge 
configuration  with  high-quality  resistors  and  connected  to  a  preamp  mounted  in  the 
can.  This  circuitry  is  used  for  controlling  the  temperature  by  regulating  the 
power  supplied  to  heating  blankets  mounted  around  each  oven. 

VSA  Processing;  System 

Six  circuit  cards  were  counted  in  a  cardcage  and  they  functioned  as  follows: 


1)  Heater  Amp  Card:  The  outputs  from  both  inner  and  outer  can  thermistors  and 


their  associated  preamps  in  the  can  are  scaled  and  amplified  on  this  card  and  then 
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are  sent  to  the  power  transistors  on  the  heat  sink  which  control  the  power  sent  to 
the  heater  blankets  around  the  ovens.  At  an  intermediate  point  in  the  signal 
routing  on  this  card,  control-voltage  signals  are  also  output  to  jacks  on  the  front 
panel  for  testing  and  monitoring  purposes.  2)  Monitor  Amp  Card:  Two  bridge/ 
amplifier  configurations  are  found  on  this  card.  These  are  connected  to  the  inner 
and  outer  monitor  thermistors.  The  outputs  are  sent  to  front-panel  janku  aid  to 
the  encoder  (via  the  analog  buffer  card).  Two  more  amplifiers  are  also  '.l.iluded 
for  measurement  of  the  resistances  of  two  other  thermistors:  one  mounted  on  the 
heat  sink,  and  the  second  on  the  mounting  plate  of  the  VSA  oven.  These  are 
intended  for  the  observation  of  ambient  temperature  within  the' gondola.  The 
output  of  this- pair  of  amplifiers  is  sent  to  the  analog  buffer  card  where,  like  the 
monitor  signals,  it  is  unaltered  and  passed  directly  to  the  encoder.  3) 

Subtract  Card:  This  board  effectively  performs  the  "subtraction"  of  the 
frequencies  of  the  two  sinusoidal  outputs  of  the  can.  One  output  is  a  square- wave 
signal  with  a  frequency  of  fd  used  in  the  Phase  Lock  Loop  (PLL)  board  in  the 
cube.  A  second  output,  another  square  wave  with  frequency  f^/N,  where  N  can 
be  selected  to  be  either  10  or  100  with  an  on-board  switch,  is  applied  to  the 
cube's  Period  Counter  (PC)  circuit.  4)  Analog  Buffer  Card:  Two  signals  from  the 
Heater  Amp  card,  the  inner  and  outer  control  voltages  are  scaled  on  f-his  board 
and  buffered  for  output  to  the  analog  section  of  the  encoder.  Likewise,  the  actual 
voltages  on  the  inner  and  outer  heater  blankets  are  scaled  and  buffered  along  with 
the  raw  battery  voltage.  Finally,  four  signals  from  the  Monitor  Amp  card,  heat 
sink  temperature  mounting  plate  temperature,  inner  monitor  temperature,  and  outer 
monitor  temperature  pass  undisturbed  through  this  card  on  their  way  to  the  analog 
inputs  of  the  encoder.  5)  Digital  Buffer  Card:  Thirty-two  digital  input  lines 
(bits)  from  the  cube  are  buffered  on  this  card,  and  30  of  them  are  output  to  the 
digital  word  inputs  of  the  encoder.  Only  the  30  Least  Significant  Bits  (LSBs)  out 
of  32  bits  are  sent  due  to  the  limitations  of  the  telemetry;  the  remaining  two  bits 
are  recovered  in  post-flight  processing.  6)  Display/CPU  Buffer  Card:  The  same 
32  digital  input  lines  as  above  are  also  input  to  this  card,  which  multiplexes  them 
in  two  separate  ways  for:  a)  the  display  device,  and  b)  the  microcomputer.  Both 
outputs  include  an  8-bit  data  word,  and  associated  address  lines.  This  card  is 
used  only  for  on-the-ground  testing  and  measurement. 

Oscillator 

The  shielded  "cube"  contains  two  circuit  boards  and  a  high-precision  10  MHz 
oscillator.  The  use  of  a  separate,  shielded  container  for  the  high-frequency  parts 
of  the  system  was  intended  to  minimize  possible  noise  problems.  The  oscillator  is 
mounted  on  a  motherboard  upon  which  is  emplaced  both  an  oscillator  oven  control 
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circuit  and  a  set  of  two  connectors  that  hold  the  PLL  and  PC  cards.  The  PLL  and 
PC  circuits  operate  concurrently.  The  purpose  of  both  is  the  high-resolution 
determination  of  f^.  The  outputs  of  both  (32  digital  signals)  are  brought  out  to 
two  Blue  Ribbon  connectors  on  the  cube.  The  pinout  of  both  connectors  is  the 
same,  so  that  the  selection  of  the  PLL  or  PC  mode  of  operation  requires  only 
connection  of  the  data  cable  (between  the  cardcage  and  the  cube)  to  the  desired 
connector.  The  PC  mode  of  operation  was  used  for  the  flight.  A  discussion  of 
the  differences  between  the  PC  and  PLL  methods  appears  below,  in  the  section  on 
individual  card  schematics. 

Encoder 

The  encoder  is  an  AFGL-provided  device  that  will  accept  three  10-bit  digital 
words  and  59  0-5  VDC  analog  signals.  These  are  multiplexed  into  one  telemetry 
stream  for  transmission. 

Heat  Sink 

The  heat  budget  of  the  gondola  is  not  a  simple  issue.  Due  to  the  lack  of 
atmosphere  at  altitude,  convection  does  not  remove  heat  from  the  gondola.  Vir¬ 
tually  all  the  heat  must  be  removed  by  radiation,  and  the  gondola  system  achieves 
equilibrium  at  a  temperature  above  that  at  1  A.  The  heat  sink  is  designed  to 
conduct  heat  out  of  the  VSA  and  move  it  to  the  gondola  frame. 

Front  Panel 

The  front  panel  is  simply  a  switch  and  connection  panel  for  ground  testing, 
and  is  completely  disconnected  during  flight. 

2. 3. 1.3  Method  for  Frequency  Counting 

Since  the  value  of  fj  at  1  g  is  roughly  64.1  Hz,  the  PLL  chip  outputs  a 
square  wave  of  256  kHz  mean  frequency.  By  its  very  nature,  the  frequency 
counting  technique  averages  the  desired  information  over  the  gate  time,  but  also 
performs  some  analog  filtering  of  the  data.  This  technique  provides  a  convenient, 
fairly  high  resolution  representation  of  the  average  frequency  of  f^,  with  the 
disadvantages  of  observation  of  frf  through  a  "fuzzy"  window  caused  by  the 
nonlinearities  in  the  analog  filtering  of  the  PLL,  and  the  phase-detector  stages  of 
the  chip.  In  order  to  provide  a  less  processed  data  set  for  the  balloon  experi¬ 
ment,  the  period  counting  method  was  actually  used. 

At  certain  times,  however,  it  was  advantageous  to  use  the  frequency  counting 
data  for  on-the-ground  testing,  since  a  direct  display  of  the  accelerations  experi¬ 
enced  by  the  sensor  (4096  *  f  )  was  available  in  real  time. 

In  period  counting,  as  in  frequency  counting,  frequency  variations  of  f^  are 
averaged  over  the  time  between  N  zero  crossings.  N  was  10  during  the  actual 
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flight,  so  that  averaging  time  at  1  g  acceleration  was  10. *(1/64.1)  or  roughly  156 
msec.  Less  filtering  would  have  been  accomplished  for  N=l,  but  then  the  data 
rate  would  have  been  too  high  for  the  acquisition  system.  With  this  method,  no 
other  filtering  was  done,  so  a  "purer"  form  of  acceleration  data  was  recorded. 

One  drawback  of  this  approach  is  the  fact  that  data  were  sampled  at  irregular 
intervals.  An  interpolation  scheme  using  bilinear  interpolation  was  used  in  order 
to  produce  a  constantly  sampled  data  stream  amenable  to  further  processing 
techniques. 

2.5. 1.4  Discussion 

Of  the  230,000  points  collected  during  the  balloon  experiment,  only  64 
"glitches"  could  be  detected  in  the  data.  Some  of  these  errors  were  attributable 
to  telemetry  dropouts,  but  some  errors  occurred  in  the  WHOI  circuitry,  in  parti¬ 
cular  the  latching  circuitry  implemented  in  the  PC  counter.  Special  care  was  taken 
in  this  circuitry  to  ensure  that  latching  would  not  occur  during  clock  transitions. 

It  would,  perhaps,  be  worthwhile  to  look  at  the  circuits  again  to  see  if  any 
improvements  can  be  made.  Nevertheless,  considering  the  environment,  an  error 
rate  of  64/250000  is  acceptable. 

The  IC  logic  type  used  throughout  was  high-speed  CMOS  (HC)  which,  with 
its  low  power  consumption  and  higher  fan  out,  seemed  to  work  very  well.  The 
only  problem  encountered  occurred  in  the  interface  between  this  type  of  logic  and 
other  more  standard  types.  For  instance,  before  the  flight,  there  was  a  flurry  of 
activity  resulting  from  a  failure  of  a  few  of  the  WHOI  chips.  These  chips  were  of 
the  HC  type,  and  were  located  at  the  interfaces  between  the  gravity  system  and 
the  encoder,  the  test  computer,  and  the  display.  The  latter  components  all  used 
the  more  standard  low  power  Schottky  TTL.  It  would  be  wise,  when  redesigning 
this  system,  to  keep  the  HC  logic,  but  to  pay  more  attention  to  interfaces  with 
other  logic  types. 

2.3.2  UNIVERSITY  OF  TORONTO  INSTITUTE  FOR  AEROSPACE  STUDIES 
(UTIAS)  INSTRUMENT  PACKAGE 

5 

The  UTIAS  flight-test  package  consisted  of  three  rate  gyros  and  three 
accelerometers.  The  rate  gyros  were  mounted  with  their  axes  orthogonal  to  each 
other,  and  would  measure  the  pendulation  and  rotation  rates  with  respect  to  a  set 
of  axes  fixed  to  the  instrument  package.  Relative  angular  displacements  would  be 


5.  DeLeeuw,  J.  H.,  and  Kung,  W.  L.  (1984)  Development  of  motion-sensing 
package  for  high-altitude  balloons,  University  of  Toronto  Institute  for 
Aerospace  Studies,  Scientific  report  submitted  under  contract 
F19628-82-C-0041 . 
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obtained  by  Integrating  the  gyro  rates.  The  vertically  mounted  accelerometer 
would  sense  the  balloon's  vertical  acceleration  superimposed  on  the  earth's  gravi¬ 
tational  field.  After  removing  the  gravity  signal  from  the  measurements,  they 
could  be  integrated  twice  to  yield  relative  vertical  displacements.  A  three-axis 
magnetometer  that  detects  the  components  of  the  earth’s  magnetic  field  about  three 
orthogonal  axes  was  also  included  in  the  balloon  motion-sensing  package. 

2. 3. 2.1  Instruments 

Rate  Gyros:  The  rate  gyros  are  GG440  A7  GNAT  miniature  rate  gyros  made 
by  Honeywell.  A  torsion  bar  spring  is  used  to  provide  a  restoring  torque  about 
the  gyro's  precessional  axis  so  that  the  precession  rate  is  proportional  to  the  input 
rate.  The  rate  of  precession  is  detected  and  converted  to  an  electrical  signal  by 
the  signal  generator /pickoff.  Two  rate  gyros,  designated  X  and  Y,  were  assigned 
to  detect  balloon  pendulation.  The  other  rate  gyro,  called  the  Z-gyro,  was  used 
to  measure  the  balloon's  vertical  rotation.  The  rate  gyros  required  400  Hz  excita¬ 
tion,  which  was  provided  by  an  Abbott  power  supply.  The  output  of  the  signal 
generator  must  be  demodulated  to  obtain  a  signal  proportional  to  input  rate. 

Accelerometers:  The  vertically  mounted  or  Z-accelerometer  should  be  able  to 
sense  the  up-and-down  motion  of  the  balloon  system.  Assuming  sinusoidal  motion, 
the  maximum  accelerations  involved  are  a  few  hundredths  of  one  g.  The  accelerom¬ 
eters  are  Sundstrand  QAllOOs.  These  instruments  are  force-feedback  acceleration 
transducers  that  feature  a  completely  elastic  seismic  suspension,  made  from  amor¬ 
phous  quartz,  that  virtually  eliminates  bearing  friction.  The  servo  electronics  are 
completely  contained  within  the  instrument,  which  also  has  a  convenient  self-test 
capability.  The  output  voltage  signal  is  developed  across  an  internal  load  resistor. 

Magnetometer:  The  three-axis  magnetometer  is  a  Schonstedt  SAM-73C.  It  is 
a  sensitive  -\strument  and  exhibits  highly  linear  response  on  all  axes.  Its  main 
purpose  is  to  provide  an  absolute  azimuth  reference  from  measurements  of  the 
horizontal  components  of  the  earth's  magnetic  field. 

2. 3. 2. 2  System  Description 

The  motion-sensing  package  hardware  is  housed  in  a  pressure-tight  aluminum 
enclosure.  Adjustable  pads  on  each  corner  and  the  two  bubble  levels  enable  the 
package  to  be  leveled  after  mounting.  The  enclosure  consists  of  a  baseplate  and  a 
cover,  with  an  O-ring  providing  the  seal  between  them.  There  is  a  port  on  the 
cover  that  allows  the  package  to  be  vacuum-tested.  A  pressure-tight  container 
was  necessary  because  the  accelerometers  were  not  hermetically  sealed.  This 
meant  that  changes  in  air  pressure  could  affect  the  instrument's  internal  damping 
and,  hence,  its  performance. 
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The  accelerometers  are  housed  in  a  triaxial  mounting  block  of  anodised  alu¬ 
minum  that  is  fastened  to  a  Lexan  base.  This  base  can  be  leveled  independently 
of  the  aluminum  baseplate.  The  arrangement  for  the  rate  gyros  is  similar.  These 
sensors  are  surrounded  on  five  sides  by  blue  styrofoam  walls  that,  with  the  Lexan 
plates,  form  an  insulated  compartment.  The  enclosure  is  heated  to  a  constant 
temperature  of  50°  C  to  prevent  external  temperature  fluctuations  from  affecting 
performance.  The  "oven"  is  heated  by  a  Darlington  power  transistor  regulated  by 
an  on /off  control  circuitry.  Two  sensors  provide  a  readout  of  the  "oven"  tem¬ 
perature,  and  the  temperature  of  the  rest  of  the  instrument  package. 

The  problem  of  maintaining  this  "oven"  at  a  constant  temperature  over  several 
hours  is  complicated  by  the  fact  that  the  rate  gyros  themselves  are  generating 
heat.  Eventually,  if  this  excess  heat  is  not  removed,  the  temperature  cannot  be 
regulated,  and  the  sensors  might  be  destroyed.  Calculations  showed  that  the 
extra  heat  could  be  absorbed  by  the  gondola  environment  if  its  temperature  re¬ 
mained  below  a  certain  level.  There  was  no  guarantee  of  this  being  the  case, 
however,  and  the  problem  was  solved  by  attaching  a  heat  sink  for  the  motion¬ 
sensing  package. 

2. 3. 2. 3  Calibration 

Magnetometers :  No  calibrations  were  performed  for  the  three-axis  magnetom¬ 
eter  because  of  the  lack  of  large  Helmholtz  coils.  Instead,  the  manufacturer’s 
calibration  data  were  used  to  reduce  the  flight  data. 

Rate  Gyros:  The  rate  gyro  calibrations  were  made  using  a  Genisco  Model 
C-181  rate  table  that  is  capable  of  rotating  at  rates  from  0  to  1200  deg/sec  in  both 
clockwise  and  counterclockwise  directions.  The  gyro  to  be  calibrated  had  its 
sensitive  axis  aligned  with  the  rotational  axis  of  the  rate  table.  This  meant  that  a 
special  mounting  jig  had  to  be  used  for  testing  the  X-  and  Y-rate  gyros.  Power 
and  signal  lines  were  routed  through  slip  rings  on  the  rate  table.  The  high-  and 
low-sensitivity  rate  gyro  outputs  were  measured  using  a  Hewlett  Packard  true  RMS 
meter.  Readings  were  obtained  as  the  rate  was  increased  in  steps  from  0  to  1 
deg/sec  and  then  decreased  for  both  clockwise  and  counterclockwise  rotation. 

Accelerometers :  Static  calibration  of  the  accelerometers  was  performed  by 
tilting  the  package  along  the  sensitive  axis  of  each  sensor,  and  measuring  the 
component  of  earth's  gravity.  The  experimental  setup  was  similar  to  that  used 
previously  at  UTIAS  for  calibrating  the  flight-test  package.  The  instrument 
package  was  fastened  to  its  mounting  jig,  and  the  whole  assembly  was  clamped  to  a 
large  turntable  that  could  be  rotated  about  a  horizontal  axis.  The  sensor  outputs 
were  read  out  on  a  true  RMS  meter. 
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3.  FLIGHT  OPERATIONS 
3.1  Preflight 

Approximately  two  weeks  ahead  of  a  scheduled  launch  date,  all  equipment  and 
personnel  are  expected  to  arrive  on  location  at  Holloman  AFB,  New  Mexico.  Equip¬ 
ment  is  unpacked  and  inspected  visually  for  damage.  Individual  modules  are 
disassembled  and  checked  for  loose  connections  or  components ,  then  reassembled 
and  tested  for  electrical  performance  following  specific  checklists.  Calibrations  are 
reverified  as  necessary  to  get  the  flight  control  integrated  with  the  ''science”  in 
preparation  for  an  "all-up"  system  checkout  run.  This  usually  consumes  about 
three  or  four  days  during  which  time  the  final  battery  preparation,  charging,  and 
packaging  are  accomplished.  After  preliminary  subsystem  tests  are  completed,  a 
full-scale,  open-loop  telemetry  and  radio  command  test  is  performed,  and  prelimi¬ 
nary  data  tapes  and  chart  records  are  run  to  verify  the  absence  of  undesirable 
interference  (or  presence  of  the  same)  that  must  be  eliminated  before  proceeding 
on  schedule. 

After  a  series  of  acceptable  preliminary  test  runs,  the  complete  payload  is 
taken  to  the  Holloman  environmental  test  facility  for  a  temperature  and  pressure- 
controlled  test  run,  lasting  from  four  to  six  hours,  to  simulat-  a  typical  ascent,, 
float,  cold  soak,  and  descent  as  would  be  encountered  during  a  balloon  flight.. 

The  chamber  controllers  attempt  to  follow  the  requested  profile,,  and  the  oppor¬ 
tunity  to  pause  and  make  changes  during  runs  permits  corrective  action  if  prob¬ 
lems  show  up  during  the  test.  This  is  a  full-scale,  open-loop,  data-taking,  test 
run.  Chamber  testing  continues  until  the  complete  system  successfully  passes  the 
full  six-hour  profile. 

After  successful  environmental  tests  are  concluded,  the  final  stages  of  flight 
preparation  begin.  This  includes  recharging  batteries,  refining  checklists  for 
prelaunch  and  flight-line  activities,  developing  a  specific  timetable,  schedule  of 
events,  and  detailed  contingency  plan  to  circumvent  potential  anomalies  should  they 
occur  during  the  launch  r  flight  phases  of  the  operation.  Once  schedules  and 
corrective  action  plans  are  finalized,  the  detailed  meteorological  analyses  of  local 
surface  weather,  flight-path  forecast,  trajectory,  and  final  impact  predictions 
begin  in  earnest,  with  daily  briefings  aimed  toward  the  commencement  of  a  final 
countdown  for  launch. 

A  typical  log  of  flight  events  (Appendices  C  and  D)  shows  a  highly  labor- 
intensive  effort  encompassing  24  consecutive  hours  of  solid  activity,  and  involving 
a  team  effort  of  between  20  to  30  people  to  produce  a  6-  to  8-hour  scientific 
balloon  mission.  The  team  effort  involves  planning,  documentation,  instrument 
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preparation,  and  testing.  It  involves  ground  support  equipment,  vehicle  readi¬ 
ness,  payload  transport,  convey  roll-out,  launch  layout,  and  flight  rigging  at  the 
launch  site,  as  well  as  final  prelaunch  verification  checkouts,  balloon  layout  and 
inflation,  launch  ascent  to  float  altitude,  change  of  float  altitude,  if  called  for,  £nd 
valve  down  and  level  off.  Following  the  completion  of  the  experiment,  it  is  nec¬ 
essary  to  set  up  and  accomplish  flight  termination. 

The  payload  descent  and  recovery  phase  can  be  critical.  When  the  scientific 
side  of  the  mission  is  finished,  the  focus  changes,  and  the  control-center  team 
concentrates  on  accurately  predicting  a  safe  and  suitable  impact  point.  They  must 
also  specify  a  location  for  descent  penetration  of  flight  level  440  (44,000  ft  MSL) 
for  transmittal  to  the  FAA  regional  control  center.  By  judicious  use  of  available 
meteorological  data,  flight-control  experience,  and  command /control  of  the  balloon 
and  payload,  the  balloon  system  is  steered  to  the  desired  termination  point  where 
the  payload  is  separated  from  the  balloon  by  radio  command.  The  balloon  is 
thereby  destroyed,  the  recovery  parachute  is  deployed,  and  a  prediction  of  para¬ 
chute  drift  to  the  impact  location  is  transmitted  to  the  tracking  aircraft  and  the 
FAA.  A  constant  plot  of  altitude  arid  position  is  maintained  throughout  flight,  and 
all  participants  are  in  constant  communication  via  several  radio  networks  and 
telephone.  The  tracking  aircraft  maintains  contact  with  the  predeployed  recovery 
mobile  force  which  is  equipped  with  suitable  vehicles,  and  tools  to  safely  handle  or 
disassemble  the  payload  if  necessary,  for  return  to  the  launch  base.  The  payload 
recovery  rate  is  exceptionally  good  and  fast.  A  refllght  is  usually  possible  after 
minor  refurbishment  and  replenishment  of  battery  power. 

3.2  Flight  History 

Flight  H83-17 

The  first  flight  of  the  gravity-measurements  experiment  was  scheduled  for  17 
September  1983  at  Holloman  AFB,  NM.  Predicted  launch  winds  were  very  favorable, 
all  prelaunoh  systems  tests  were  "go,"  and  the  balloon  was  inflated  at  Nenninger 
Site  on  Holloman  AFB,  NM.  Immediately  after  the  balloon  bubble  was  released, 
however,  the  launch-arm  operator  observed  a  long  tear  in  the  uninflated  material, 
and  the  balloon  developed  a  large  sail.  The  crane  holding  the  gondola  backed 
away  from  the  launch  arm  and  dropped  the  gondola.  The  sailing  balloon's  lift  in 
one  direction,  and  the  opposite  pull  of  the  crane,  had  resulted  in  very  high 
tension  at  the  hydraulic  boom  holding  the  gondola.  Subsequent  investigation 
showed  that  under  that  tension  the  boom  would  extend  several  feet,  causing  the 
payload-release  mechanism  to  open.  Although  the  shell  of  the  cylindrical  gondola 
was  badly  damaged,  it  had  provided  excellent  protection  for  the  instrumentation 
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inside;  all  of  the  sensors  and  electronics  survived  unimpaired.  The  crane  mal¬ 
function  was  corrected,  the  gondola  was  refurbished,  and  the  flight  was  resched¬ 
uled  for  October  1983. 

Flight  H83-19 

All  of  the  instrumentation  was  thoroughly  retested  and  reassembled  in  the 
gondola.  Early  in  the  morning  of  11  October  1983  the  preflight  checkout  was 
completed  smoothly,  and  the  ballopn  and  gondola  (Figures  3.1  and  3.2)  were 
transported  to  the  Nenninger  (launch)  Site.  As  the  0600  Local  Time  launch  hour 
approached,  the  Pibal  runs  at  15-min  intervals  were  indicating  that  although 
surface  winds  were  not  exceeding  the  10  to  12  knot  limit  for  this  launch,  the 
low-level  winds  were  increasing  to  unacceptable  speeds.  The  flight  meteorologist 
predicted  a  short  interval  of  relative  calm  shortly  after  sunrise,  and  inflation  was 
delayed  accordingly.  The  predicted  improvement  in  low-level  winds  did,  in  fact, 
occur,  the  balloon  was  inflated  (Figure  3.3),  and  the  launch  maneuver  was  skill¬ 
fully  executed  at  0811  Local  Time,  with  the  crane  running  the  full  extent  off' the 
paved  site.  The  ascending  balloon  (Figure  3.4)  traveled  approximately  southeast 
until,  just  before  it  entered  float  phase  at  1004  Local  Time,  it  abruptly  changed 
heading  and  continued  roughly  northeast  until  the  flight  was  terminated.  The 
flight  path  is  shown  in  Figure  3.5.  Two  hours  of  data  were  recorded  at  the 
100,000  ft  float  level;  then,  by  commanded  gas  valving;  and  some  judicious  pouring 
of  ballast,  the  balloon  was  nudged  into  the  83,000  ft  level.  Its  float  motion  there, 
however,  was  somewhat  disturbed,  and  at  1430  Local  Time  the  termination  command 
was  given  to  separate  the  gondola-parachute  train  from  the  balloon.  The  para¬ 
chute  recovjry  system  had  purposely  been  designed  for  low-g  impact,  using  a 
100— ft  diam,  flat  circular  parachute;  consequently,  the  descent  time  was  unusually 
long.  The  descent  was  observed  by  radar  down  to  35,000  ft  at  2015.  Ground 
impact  was  recorded  by  the  chase  aircraft  at  2124.  The  gondola  landed  in  a 
plowed  field,  just  off  a  road,  southeast  of  Lovington,  NM,  and  was  recovered  in 
excellent  condition. 


4.  DATA  ANALYSIS 
4.1  Objectives 

The  primary  objective  of  a  feasibility  experiment  is  to  answer  the  question 
"Can  it  be  done?"  The  other  objectives  generally  include  an  attempt  to  actually 
"do  it"  as  well  as  possible,  with  appropriate  caveats  and  limitations.  In  this 
project,  a  great  deal  of  data  were  collected,  organized,  inspected,  and  analyzed  by 
several  different  people  and/or  groups.  Since  this  first  flight  was  a  feasibility 
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Figure  3.1  A  Photograph  of  the  Payload,  named  DUCKY  la,  Just  After  Transport  to  the 
Launch  Site 
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flight,  the  results  are  neither  rigorous  nor  conclusive.  With  this  view,  assembling 
the  existing  data  analysis  into  a  coherent  and  organized  package  is  difficult.  The 
current  objectives  fall  into  two  broad  groups;  1)  How  good  are  the  data,  and 
what  should  be  done  on  subsequent  flights  to  improve  the  data  quality?  and  2) 

What  is  the  measured  gravity,  and  how  do  the  measurements  compare  with  mathe¬ 
matical  models? 

An  overview  of  the  data  is  shown  in  Figures  4.1  and  4.2.  The  former  shows 
the  altitude  profile  and  raw  gravimeter  data  for  the  entire  flight;  the  latter  shows 
the  same,  but  only  during  float.  Note  that  there  are  two  "quiet"  periods,  just 
before  11  A.M.  and  just  before  12  P.M.,  local  time.  Most  of  the  data  analysis  that 
follows  focuses  on  the  second  quiet  time  as  it  also  showed  the  least  rotation  about 
the  vertical  axis.  The  strategy  is  to  understand  a  short  segment  of  "quiet"  data 
before  attempting  to  analyze  the  entire  flight.  In  the  end,  enough  was  learned 
from  this  analysis  to  justify  and  contribute  to  a  second  flight. 

The  most  important  point  for  understanding  the  data  analysis  issues  is  that 
this  really  is  a  complex  navigation  problem.  A  complete  solution  requires  a  thor¬ 
ough  understanding  of  not  only  the  geophysics  of  gravity  and  gravity  modeling, 
but  also  of  navigation  systems  and  analysis  techniques.  The  analysis  to  be  pre¬ 
sented  here  is  intended  to  illuminate  the  obvious  problems  and  solutions  to  be 
addressed,  and  does  not  come  anywhere  near  a  complete  solution.  Part  of  the 
reason  for  that  is  due  to  the  great  amount  of  effort  that  was  required  just  to 
assess  the  raw-data  quality  and  assemble  that  data  into  an  organized  investigator 
tape;  and  part  is  due  to  the  severe  limitations  resulting  from  unsatisfactory 
tracking  data  that  undermined  the  only  "ground  truth"  data  available.  Neverthe¬ 
less,  significant  progress  was  made  resulting  in  improved  understanding  of  balloon 
dynamics,  gravity  measurements  (although  not  to  required  accuracy),  and  the 
improved  design  of  the  gondola  and  instruments  for  subsequent  flights. 

The  data  analysis  divides  into  roughly  two  parts;  corrections  for  translational 
motions,  and  corrections  for  rotational  motions.  The  translational  motions  are 
derived  from  the  radar-tracking  data,  and  have  been  "smoothed"  in  the  vertical 
direction  using  a  Kalman  Filter  with  the  vertical  accelerometer  data  as  additional 
input.  This  procedure  used  the  higher  frequency  components  of  the  accelerometer 
data  only,  minimizing  the  low-frequency  components.  Thus,  we  believe  we  have 
avoided  the  "chicken-and-the-egg"  problem,  since  it  is  the  lowest  frequency  com¬ 
ponent  of  vertical  acceleration  that  we  are  after.  Rotational  corrections  were  not 
made  due  to  lack  of  appropriate  navigation  software.  However,  the  gyro  and 
magnetometer  data  were  inspected  to  give  an  indication  of  vertical  angle  errors, 
and  thus  the  stability  of  the  gondola.  The  direct  purpose  of  such  an  inspection  is 
to  determine  the  need  .for  a  stabilized  platform  in  future  flights. 
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Comparison  of  Flight  Altitude  and  Gravimeter  Measurements  for  the  Entire  Flight 


Figure  4.2  A  Comparison  of  Flight  Altitude  and  Gravimeter  Measurements  at  Float  Altitude  Only 


4.2  Prediction  Filter  to  Fill  in  Gaps  of  Radar-Tracking  Data 


4.2.1  INTRODUCTION 

One  of  the  problems  encountered  during  the  post-flight  data  analysis  was  the 
existence  of  numerous  data  gaps  in  the  WSMR  radar-tracking  data.  The  presence 
of  gaps  in  the  data  can  cause  severe  problems  for  tracking  algorithms ,  and  may 
have  caused  problems  for  DMA  in  deriving  gravity  estimate  values  from  tracking 
data.  For  these  reasons,  an  algorithm  was  developed  to  predict  the  missing  data 
to  preserve  its  continuity. 

Sixty-one  gaps  were  identified,  ranging  from  2  to  30  samples  in-length.  In 
addition,  three  large  gaps  were  identified  that  were  259,  199,  and  839  samples 
long.  Because  of  their  length,  these  gaps  were  not  filled  since  the  performance  of 
any  algorithm  would  be  highly  suspect  over  such  long  periods. 

In  the  next  section  we  present  the  algorithm  used  to  predict  the  missing 
data.  This  algorithm  is  based  on  the  theory  of  linear  prediction.  Following  a 
discussion  of  the  algorithm,  a  brief  discussion  of  its  implementation  and  some 
graphical  results  are  presented. 

4.2.2  LINEAR  PREDICTION  ALGORITHM 

The  approach  taken  to  fill  in  the  gaps  is  to  use  the  technique  of  linear 
prediction.  This  technique  assumes  an  auto-regressive  (AR)  time  series  model  for 
the  data,  and  uses  this  model  to  predict  the  missing  data  based  on  the  existing 
data  surrounding  the  gap.  The  signal  x(n)  is  assumed  to  be  a  linear  combination 
of  past  values  and  some  input  u(n): 

P 

x(n)  =  E  a-x(n-k)  +  u(n)  (4.1) 

k=l  K 


2 

The  input  u(n)  is  assumed  to  be  a  white  noise  term  with  variance  a  .  The  model 
order  P  must  be  estimated,  and  the  AR  coefficients  {a^}  must  be  determined. 
Experiments  with  the  tracking  data  showed  that  P=8  is  a  good  model  order  to  use. 

The  sample  values  x(nQ),  x(no+l),...,x(n1)  are  assumed  to  be  missing. 
Thus,  the  gap  is  L=nj-n  +?  points  wide. 

The  data  immediately  preceding  the  gap  are  used  to  estimate  the  forward 
prediction  coefficients  a^  via  the  Burg  method.  This  method  is  part  of  the  basic 

6.  Makhoul,  John  (1975)  Linear  prediction:  a  tutorial  review,  Proc.  IEEE  63 
(No.  4). 

i.  Burg,  John  P.  (1975)  Maximum  Entropy  Spectral  Analysis.  Ph.D.  Thesis, 
Stanford  University. 
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I*S*P  package  under  the  function  name  LNP.  Five-hundred  points  are  used  to 
form  the  estimates.  These  prediction  coefficients  are  used  in  the  following  formula 
to  yield  the  forward  prediction  data: 


A 

Xf 


(n) 


I  a.x,(n-k), 
k=l  K  1 


n=n 


o’ 


V1 


(4.2) 


where  x^(n)=x(n)  for  n<nQ. 

Similarly,  the  data  immediately  following  the  gap  are  used  to  estimate  the 
backward  prediction  coefficients  {b^}.  Again,  the  Burg  method  is  used  and  500 
points  are  used  to  form  the  estimate.  The  backward  prediction  coefficients  are 
used  in  the  following  formula  to  form  the  backward  prediction  data: 


A  /  X 

xb(n) 


?  b,_x 


k=l 


kAb(n+k), 


n=n^,  n^-1,  ...,  nQ 


(4.3) 


where  xb(n)=x(n)  for  n<nj. 

The  final  prediction  estimate  x(n)  is  formed  by  adding  weighted  versions  of 
the  forward  and  backward  estimates.  The  weights  are  proportional  to  the  error  in 
the  predictors,  where  the  error  will  increase  as  n  is  farther  **way  from  nQ  for  the 
forward  predictor,  and  n^  for  the  backward  predictor.  The  final  prediction 
equation  is  then  given  by: 


x(nQ  +  m)  = 


L-l-m 

Y  ( n  +m  ^  + 

m 

L-l 

L-l 

xb(n0+m). 
m=0 ,  1,  ...,  L— 1. 


(4.4) 


4.2.3  IMPLEMENTATION  ISSUES 

If  the  data  have  a  DC  offset  and/or  linear  trend,  they  can  cause  errors  in 
estimating  the  prediction  coefficients,  and  will  affect  the  prediction  results  ad¬ 
versely,  Therefore,  the  DC  offset  and  trend  are  removed  from  the  data  preceding 
the  gap  and  from  the  data  following  the  gap  before  estimating  the  prediction 
coefficients,  and  then  forward  and  backward  prediction  coefficients  are  computed 
from  the  detrended  data.  Detrended  forward  and  backward  predictions  are  com¬ 
puted  using  Eqs.  (4.2)  and  (4.3).  Then  the  DC  offset  and  trend  are  added  back 
to  the  predictions,  and  the  weighting  in  Eq.  (4.4)  is  finally  applied. 


31 


V  V  V  V 


One  problem  that  was  encountered  during  the  graphical  inspection  of  the  gaps 
was  that,  very  often,  the  data  immediately  surrounding  the  gaps  seemed  to  be 
suspect.  This  is  shown  in  Figures  4.3,  4.4,  and  4.5  where  the  gap  is  indicated 
by  the  two  vertical  lines.  Note  that  Figures  4.4  and  4.5  show  two  gaps  in  the 
time  span  shown.  These  figures  show  that  the  data  just  before  the  gap  and  just 
after  the  gap  are  noisier  than  data  further  away  from  the  gap.  For  this  reason  it 
was  decided  to  extend  the  size  of  the  gap  (or  to  predict  more  points)  for  gaps 
that  showed  this  tendency.  Also  displayed  in  Figures  4.3,  4.4,  and  4.5  are  the 
corresponding  predictions.  These  are  the  smooth  curves  that  pass  over  the  gaps. 

In  Figures  4.4  and  4.5  where  two  gaps  are  so  closely  located,  the  two  are  treated 
as  a  single  gap. 

4.3  Kalman  Filter  Model  for  Smoothed  Eatinutea  of  Balloon  Position 

In  Section  4.2,  we  considered  a  Kalman  Filter  (KF)  model  for  estimating 
balloon  velocity  and  acceleration  using  radar-tracking  data  only.  Inputs  to  this 
model  were  the  radar  position  tracking  observations  and  associated  errors,  which 
gave  rather  unbelievable  results  for  velocity  and  acceleration.  In  looking  at  some 
of  the  observations  we  concluded  that  there  was  little  chance  that  the  balloon  was 
actually  doing  what  the  observations  implied,  and  that  tracking  position  errors 
given  to  us  on  tape  were  overly  optimistic.  The  reported  errors  were  on  the 
order  of  2  to  4  m,  but  we  felt  they  were' much  higher.  Using  these  small  tracking 
errors  in  the  KF  model  meant  that  we  assumed  the  observations  were  quite  good, 
leading  to  less-than-perfect  velocity  and  acceleration  results.  In  this  problem  we 
use  the  KF  again,  but  this  time  a  heuristic  estimate  is  made  of  the  position  errors, 
given  the  data,  and  the  high-gain  Z-accelerometer  observations  are  incorporated  as 
well.  The  attempt  here  is  to  reprocess  the  data  and  look  at  the  smoothed  position, 
velocity,  and  acceleration  estimates,  given  observations  of  both  balloon  position 
and  acceleration. 

4.3.1  KALMAN  FILTER  MODEL 

Define  the  state  vector 

x(k)  =  [x(k)  *(k)  Jt(k)  ]T  (4.5) 

where  x(k)  is  the  balloon  position  at  time  k,  and  £(k)  and  K(k)  are  the  first  and 
second  derivatives,  respectively.  The  state  model  for  the  KF  is  based  on  constant- 
acceleration  model  with  random  perturbations: 
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Figure  4.3  Example  of  a  Gap  of  Eight  Points  in  WSMR  Tracking  data.  Gap  is  indicated 
by  the  two  vertical  lines  in  the  plot.  Also  shown  is  the  predicted  data,  which  is  the 
smooth  line  through  the  gap.  Note  that  the  predicted  data  actually  spans  more  points 
than  the  gap  alone ,  due  to  the  larger  fluctuations  in  the  data  surrounding  the  gap 


27940 
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Figure  4.4  Example  of  Section  of  Data  With  Two  Gaps,  each  of  Length  6.  Each  gap  is 
indicated  by  two  vertical  lines  in  the  plot.  Also  shown  is  the  predicted  data,  which  is 
smooth  line  through  the  gaps.  Note  that  the  predicted  data  actually  spans  more  points 
than  the  gaps  alone,  due  to  the  larger  fluctuations  in  the  data  surrounding  the  gaps 


(sj»j»in)  saniinv 
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Figure  4.5  Example  of  Section  of  Data  With  Two  Gaps,  Each  of  Length  6.  Each  gap  is 
indicated  by  two  vertical  lines  in  the  plot.  Also  shown  is  the  predicted  data,  which  is 
the  smooth  line  through  the  gaps.  Note  that  the  predicted  data  actually  spans  more  points 
than  the  gaps  alone,  due  to  the  larger  fluctuations  in  the  data  surrounding  the  gaps 


x(k+l)  =  Fx(k)  +  Gu(k) 


(4.6) 


0 

i 

.  1 

A  .5A2| 

where  G  = 

0 

[and  F  = 

1  0 

1  A  ! 

1 

* 

1 

1° 

0  a 

(system  matrix)  (4.7) 


where  A  is  the  sampling  interval,  a  is  the  acceleration  parameter  defining  the 
bandwidth  of  acceleration,  and  u(k)  is  the  driving  white  noise. 


q  =  Q  *  E[u2(k)] 


(4.8) 


is  the  white  noise  power.  We  are  dealing  with  A  =  .1  sec  and  the  parameters  a 
and  q  will  be  estimated  from  the  data. 


Observation  Equation: 

We  observe  the  balloon  position  x(k)  and  acceleration  X(k)  so  these  are  our 
observations  z(k). 


z(k)  = 


i  x(k) , 

n  ....  !j 

!  I  X(k)  I 


=  Hx(k)  +  v(k) 


(4.9) 


and 


H  = 


10  0 
0  0  1  ! 


(4.10) 


where  v(k)  is  observation  noise.  Define 


R(k)  =  E[v(k)vT(k)J  =] 

I 


1 

r2i 


(4.11) 


In  our  first  KF  attempt  we  used  r^=4  (from  the  tracking-error  tape),  but  this  time 
we  will  make  a  heuristic  estimate  of  and  r2. 


The  Kalman  Filter: 

The  KF  performs  a  five-step  iteration  for  each  time  k: 

(1)  State  estimate  extrapolation: 

\xo(k+l)=Fx(k)  (4.12) 
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where  tx<k)  -  state  vector  estimate 

\x(ki-i)  -  slate  vector  estimate  extrapolation 

(2)  Error  covariance  extrapolation: 

P0(k+l)=FP(k)FT+GQGT  (4,13) 

where  P(k)=E[  (x(k)-x(k)(x(k)-x(k))T]  (4.14) 

is  the  error  covariance  matrix,  and  PQ(k+l)  is  the  error  covariance  matrix 
extrapolation 

(3)  Kalman  gain  computation: 


K  =  Po(k+l)HTtHP0(k+l)HT+R]“T 


where  K  is  the  Kalman  gain  matrix, 

(4)  State  estimate  update: 

x(k+l)=x^(k+l)+K  (z(k+l)-Hjx\(k+l)  ] 

(5)  Error  covariance  update: 


(4.15) 


(4.16) 


P(k+l)=[I-KH]Po(k+l).  (4.17) 

The  block  diagram  form  is  shown  in  Figure  4.6. 

4.3.2  DATA  USED  AND  PARAMETER  SELECTION 

The  data  used  were  a  section  of  data  from  1730-1807  zulu.  The  original  20 
Hz  Z-tracking  data  were  decimated  to  10  Hz  for  this  problem.  Plots  of  the  de¬ 
trended  and  DC  bias-removed  observations  are  shown  in  Figure  4.7.  Note  that  in 
the  lower  figure  the  detrended  Z-accelerometer  data  were  multiplied  by  -1  to  make 
tiie  directions  of  the  Z-axes  the  same.  To  select  a  and  g,  we  tried  a  few 
parameters,  synthesized  some  data,  and  settled  on  the  values: 


a  =  0.999  - 

q  =  2.5  *  10-'  (4.18) 


To  select  r^  and  r2  we  looked  at  the  observations.  We  selected  a  value  of  r^  = 
225  =  152  m2  for  the  tracking  data,  and  r.,  =  5xl0-7  =  (.0007  m/sec2)2. 
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4.3.3  RESULTS 

The  results  of  the  XF  output  are  shown  in  Figures  4.8,  4.9,  and  4.10.  In 
Figure  4.8  the  tracking  observations,  x(k),  are  plotted  against  the  position  esti¬ 
mate  x(k)  for  the  first  204.8  sec  of  the  file.  Figure  4.9  X(k)\  and, ^  X(k)\ are 
plotted  for  the  first  204.8  sec  of  the  file.  Finally,  in  Figure  4.10,  'X(k)  is  plotted 
for  the  first  204. C  sec. 

4.4  Recursive  Estimation  Approach  to  Signal  Analysis 

As  part  of  the  exploratory  development  effort,  the  balloon  data  were  sub¬ 
jected  to  a  brief  analysis  using  the  traditional  approach  of  Wiener  Filtering,  or 
least-mean-square  estimation.  For  a  very  large  number  of  data  points  examined 
sequentially  in  time,  the  Wiener  Fitter  can  be  constructed  as  a  recursive  digital 
algorithm.  In  the  nonrecursive  case,  the  Wiener  Filter  problem  can  be  stated  as 
follows: 

Given  a  set  of  data  points  x^,  j=l,  2,  ...,  m,  where  each  Xj  consists  of  a 
signal  S  (in  this  case  constant  in  time)  plus  additive  uncorrelated  zero-mean  noise 
n^,  we  wish  to  obtain  a  "best"  estimate  of  the  signal,  S,  which  can  be  written  in 
general  form  as 

A  m 

S  =  £  hx  (4.19) 

j=l  1  * 

A 

in  other  words,  S  is  a  linear  combination  of  all  available  data,  and  the  problem  at 
hand  is  to  construct  an  algorithm  that  will  generate  the  weighting  coefficients  hj 
such  that  the  S  so  obtained  is  the  "best"  estimate  of  the  signal  S,  with  "best"  left 
undefined  in  the  most  general  sense.  In  the  Wiener  Filter  approach,  "best"  is 
defined  as  that  estimate  of  S  which  minimizes  the  mean-square-value  of  the  error  e: 

S[e2]  =  E[(S  -  S)2]  =E[(S  -^Yx.)2]  |  (4.20) 

j=l 

The  usual  approach  is  to  differentiate  the  above  expression  with  respect  to  each  hj 
and  set  equal  to  zero,  resulting  in  a  set  of  m  equations  in  the  unknowns  h^: 

E[ex.]  =0;  j  -  1,  2,  ...,  m.  (4.21) 
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Figure  4.8  Comparison  of  Raw  Radar  Altitude  Data  With  the  Kalman  Filter  Output 


0.015 


'ignre 


fith  the  Kalman  Filter  Estimate 


The  resulting  matrix  equation  would  have  to  be  solved  to  obtain  the  coefficients 
2 

that  minimize  E[e  J.  The  set  of  coefficients  becomes  the  Wiener  Filter  for  this 

problem,  and  we  can  say  that  the  data  has  been  ''passed  through  a  Wiener  Filter" 

A 

to  recover  the  best,  or  optimal,  estimate  S. 

For  a  large  data  set  such  as  the  balloon  experiment,  this  approach  would  be 
unwieldy  in  the  extreme.  A  variation  of  this  approach,  however,  lends  itself  to 
optimal  estimation  using  many  noise-corrupted  data  points  in  sequence,  that  of 
recursive  least-squares  estimation.  In  this  approach,  the  estimate  is  updated  as 
each  new  data  point  is  obtained,  with  the  new  estimate  being  the  old  estimate  plus 
a  weighted  fraction  of  the  difference  between  the  old  estimate  and  the  new  data: 

\+i  =  sk +  Wxk+1 (4-22) 

A 

where  the  obvious  problem  now  is  to  obtain  values  for  each  bv  such  that  Sk  is  the 
8  K 

optimal  estimate. 

Deriving  the  functional  form  of  b^  is  dependent  c  '  the  characteristics  of  the 
physical  process  being  modeled.  In  particular,  if  S  is  a  time-independent  random 
variable,  that  is,  its  value  during  data  acquisition  is  constant  but  unknown,  then 
recursive  Wiener  Filtering  can  be  used,  provided  an  appropriate  form  for  b^  can 
be  found.  If,  on  the  other  hand,  the  signal  is  time-dependent  (written  S(tk)), 
then  the  optimal  estimation  process  is  embodied  in  the  Kalman  Filter,  which  can  be 
written: 

St+1  =  *t+1,t  St  +  GtXt>  (4.23) 

which  corresponds  to  the  time-invariant  recursive  equation  above  with  the 
substitutions 


*t+l,  t  =  1  -  bk+1  and  Gt  =  bR+1  (4.24) 

In  the  Kalman  Filter  notation,  S  and  x  are  system  state  and  input  vectors,  Gt  is 
the  input  matrix  (set  of  input  gain  factors),  and  *t+1  t  is  the  transition,  or 
system,  matrix  describing  how  the  system  state  has  changed  from  data  point.  x(k) 
to  data  point  x(k+l).  Gt  and  t  are  written  with  time  subscripts  to  indicate 
that,  unlike  b^,  they  are  not  time-invariant. 

8.  Schwartz,  M.,  and  Shaw,  L.  (1978)  Signal  Processing,  McGraw-Hill,  New 
York,  Ch.  7. 
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In  the  absence  of  an  analytic  physical  model,  the  form  of  the  recursive 

o 

weighting  factor  is  notably  simple.  The  derivation  for  b^  will  not  be  given 
here;  the  functional  form  for  a  zero-slope  estimate  is  simply 


k  k  +  b  k  +  ijJ/EtS2]  (4.25) 

2  2 
where  o  is  the  variance,  or  power,  of  the  additive  noise  and  E[S  ]  is  the 

n  2 
variance  of  the  signal.  For  the  gravimeter  data  application,  howevers,  E[S  1  is 

2 

essentially  zero,  while  on  is  quite  large  and  measurable.  Signal  power  is,  of 

course,  close  to  zero,  because  the  true  value  of  g  changes  very  slowly,  and  in 

very  small  fractional  amounts,  as  the  balloon  describes  its  trajectory.  Therefore, 

2  2 

the  noise-to-signal  ratio  b  =  on/E[S  ]  is  extremely  large,  and  a  rough,  practical 
approximation  must  be  used  for  recursive  estimation  to  be  of  any  utility  in  this 
application. 

In  exploring  possible  processing  approaches,  various  forms  of  b^  were 
evaluated,  but  none  were  found  to  be  especially  useful  for  the  problem  at  hand. 
Specifically,  the  denominator  of  b^  was  adjusted  by;  (1)  assigning  a  heuristic 
value  to  b,  the  noise-to-signal  ratio;  (2)  including  a  term  to  reflect  signal  slope  in 
the  neighborhood  of  x^;  (3)  including  an  additional  term  to  reflect  the  deviation  of 
xk  from  the  signal  mean,  normalized  by  ox;  and  (4)  several  combinations  of  the 
three.  Due  largely  to  the  nature  of  the  data,  however,  (the  signal  being  a  con¬ 
stant  "bias"  with  additive  noise),  recursive  estimation  gave  results  similar  to 
simply  calculating  the  mean  of  the  data.  Each  v.ata  point  used  for  both  mean  and 
recursive  estimation  gave  results  similar  to  simply  calculating  the  mean.  Each  data 
point  was  corrected  previously  for  two  measurable  error  sources,  altitude  and 
Eotvos  effect.  The  data  was  sampled  at  a  10  Hz  rate,  forming  10  files  of  6144 
samples  (614.4  sec)  on  each  of  four  channels:  (1)  VSA  output,  (2)  altitude, 

(3)  calculated  Eotvos  effect,  and  (4)  reference  DMA  gravity  estimates  for  the 
position  at  that  time.  The  VSA  and  DMA  reference  channels  were  adjusted  for  a 
reference  altitude  of  30.300  km,  with  Eotvos  applied  to  the  VSA  channel  only. 

Thus,  the  total  span  of  the  data  was  102.4  min  in  length;  the  recorded  telemetry 
from  1630Z  to  1810Z. 

The  results  of  this  preliminary  processing  effort  are  tabulated  in  Table  4.1. 
Means  and  variances  for  each  of  the  four  channels  are  shown,  plus  recursive 
estimation  results  for  three  different  forms  of  the  coefficient  b^.  The  main  points 
of  this  analysis  are. 

(1)  Both  the  average  and  recursive  estimate  for  individual  files  of  10.24  min 
in  length  show  a  very  high  level  of  variability;  the  average  value  of  gravity  g 
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Table  4.1.  Results  From  10  Consecutive  10-sec  Segments.  Altitude  is  derived  by 
Kalman  smoothing  of  radar  data  and  flight  accelerometers.  Eotvos  Is  computed  from 
radar  data  only.  DMA  gravity  models  are  shown  on  raw  form  (g,),  then  continued 
to  30.3  km  altitude  and  970,000  mGals  subtracted  (g,).  VSA  data  are  shown  on 
raw  form  (VSA),  with  970,000  mGals  subtracted  (VSA„)  and  then  continued  to  30.3 
km  altitude  with  Eotvos  corrections  applied  (VSA,).  Tn  principle,  g,  and  VSA, 
show  data  as  though  the  balloon  were  stationary  at  30.3  km  altitude  with  a  970;000 
mGals  bias  removed.  Instrument  and  processing  noise  is  clearly  evident.  Although 
there  are  a  number  of  ways  to  interpret  these  results,  the  best  is  to  look  at  the 
100-sec  average 


File 

Altitude 

Eotvos 

g.  (DMA) 
raw 

g,(DMA) 
tg-970), 
30.3  km 

VSA 

raw 

VSA, 

(g-970; 

VSA, 
(g— 9707 , 
30. 3  km, 
Eotvos 

m 

mGals 

mGals 

mGals 

mGals 

mGals 

mGals 

<7 

cr 

cr 

1 

0 

a 

1 

30415 
— 57 

-43.1 

“50 

970230 
- 45 

226.5 

3.9 

970142 

2553 

136.6 

215.4 

TO 

2 

30506 

6? 

-57.3 

“50 

970231 
- 55 

230.6 

4.0 

970000 
— 555 

3.0 

123.7 

571 

3 

30523 

-82.8 

970212 

226.5 

970144 

151.6 

303.3 

14 

54 

3.5 

455 

477 

4 

30490 

-84.6 

970218 

229.5 

969980 

-24.9 

118.3 

2? 

TO 

48 

4.1 

535 

575 

5 

30400 
— U 

-63.2 

TO 

970226 
- 34 

227.9 

4.2 

970129 
— 5TI 

128.2 

222.2 

“O 

6 

30288 
— 35 

-70.2 

TO 

970285 

58 

223.9 

“O 

920113 

975 

107.9 

176.4 

“ O 

7 

30256 
— 35 

-73.8 

TO 

970290 

47 

230.1 

or 

969999 

1375 

-14.0 

46.3 

TO 

8 

30177 

-73.8 

970304 

221.3 

970328 

332.8 

368.5 

23 

T 777 

46 

~T75 

1026 

TO 

9 

30235 

-77.6 

970291 

225.3 

970225 

232.1 

290.0 

23 

TO 

47 

~4T 

1763 

1.7 

10 

30208 

78 

-92.4 

TO 

970286 

36 

234.8 

7T 

969973 
’ — 554 

-30.8 

33.3 

6.8 

Sum 

227.6 

189.7 

~T72 

8.0 

Average  over  10  segments  (100  seconds  total)  [g(DMA)-VSA]  at  30.3  km  =  37.  9  + 
6.3  mGals 
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(with  a  bias  of  970  Gals  removed)  ranges  from  a  low  of  0.3325  Gals  to  a  high  of 
0.36852  Gals,  the  gross  average  being  0.1897  Gals.  However,  when  all  10  files  are 
linked  and  recursive  estimation  is  used  for  the  entire  data  stream  of  102.4  min, 
the  variability  drops  to  extremes  of  0.1945  Gals  to  0.2842  Gals,  with  a  final  esti¬ 
mate  n?  0.1963  Gals  for  o ovsmple. 

(2)  The  best  results  with  recursive  estimation,  shown  in  (1)  above,  were 
obtained  by  including  an  additional  b-factor  in  the  correction  coefficient  denomi¬ 
nator,  being  the  product  of  the  average  absolute  value  of  local  slope  times  one, 
plus  the  deviation  of  that  data  point  from  the  segment  mean,  both  quantities 
normalized  by  the  segment  standard  deviation,  the  product  then  assigned  an. 
exponent  of  8.  This  additional  factor  appeared  to  attenuate  the  noisiest  segments 
of  the  data  stream  sufficiently  to  allow  a  small  improvement  in  the  final  result 

for  'g'. 

(3)  The  theoretical  value  of  'g'  was  obtained  two  different  ways  and  was 
found  to  be  in  very  good  agreement.  The  average  of  the  local  height  corrected 
DMA  values,  which  include  the  upward  continuation  of  surface  anomalies,  and  a 

Q 

calculation  of  average  'g'  at  altitude  using  Table  11  of  OSU  Report  #360,  both 
yielded  a  value  of  0.227  for  'g'  (970.227  Gals).\  This  placed  our  preliminary 
experimental  results  approximately  30  to  35  mGals  low  in  comparison  with  theory. 
This  error  could  easily  be  accounted  for  by  a  slight  vertical  misalignment  of  the 
VSA,  or  by  a  bias  in  the  electronics. 

(4)  An  examination  of  Table  4.1  shows  that  the  variance,  or  noise  power,  can 
be  assigned  almost  entirely  to  the  VSA  channel,  which  on  average  is  about  600 
times  the  geometric  average  of  the  variances  of  the  altitude  and  Eotvos  channels, 
the  two  t  - .  g  comparable  in  noise-power  content.  An  autocorrelation  of  the  VSA 
channel  reveals  that  virtually  all  the  noise  power  lies  in  a  narrow  frequency  band 
centered  about  a  280-sec  oscillatory  period,  very  close  to  the  characteristic  period 
of  vertical  oscillations  of  a  fully  inflated  balloon  in  an  isothermal  atmosphere, 
typically  calculated  at  about  255  sec.*®  Therefore,  improved  signal  conditioning 
emphasizing  narrowband  noise  smoothing  should  be  of  primary  consideration  in 
future  flight-data  analysis. 
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’  4.j  Rotation 


4.5.1  INTRODUCTION 

All  formal  data  analysis  described  thus  far  dealt  with  translational  motions 
only,  and  are  a  summary  of  the  best  of  many  approaches  tried.  A  complete  solu¬ 
tion  requires  proper  application  of  strap  down  inertial  guidance  concepts  and 
gravity-data  analysis.  Since  rotational  analysis  is  complex,  and  mutual  interaction 
among  the  different  data  sets,  (including  translational  analysis),  is  nonlinear,  such 
a  complete  solution  does  not  exist  yet.  Short  of  processing  all  the  data  simulta¬ 
neously,  such  as  with  an  all-inclusive  Kalman  Filter,  a  rational  approach  requires 
solving  the  translational  problems  first,  then  approaching  the  rotational  problems. 
The  latter  is  described  here  qualitatively.  Specifically,  such  insight  is  needed  to 
design  the  data  analysis  scheme  properly  for  subsequent  flights. 

4.5.2  INSTRUMENTS  AND  METHODS 

Three  sets  of  motion-sensing  instruments  were  flown  in  addition  to  the  gravi¬ 
meter.  They  arc  the  accelerometers,  the  rate  gyros,  and  the  magnetometer;  all 
three-axis  instruments.  All  give  useful  information  about  gondola  rotational  motions. 
As  high-frequency  (greater  than  0.1  Hz)  noise  has  not  yet  been  properly  character¬ 
ized,  a  0.1  Hz,  Low-Pass,  Zero-Phase  Shift  Butterworth  Filter  was  applied  to  the 
data.  Thus,  the  principal  signals  evident  are  of  periods  longer  than  10  sec.  Raw 
data  will  be  shown  in  the  figures  for  background  information  only. 

Accelerometers 

The  three-axis  accelerometer  on  board  is  oriented  with  Z  up,  X  out  the  front 
door,  and  Y  the  orthogonal  axis  in  a  right-handed  system.  We  start  with  the 
concept  of  the  "true  vertical,"  which  is  the  local  perpendicular  to  the  geoid  and  is 
the  locally  stable  direction  which  the  gondola  system  sees  as  "down."  Since  the 
gondola  system  has  alignment  errors  (static),  and  kinematic  motions  (dynamic),  the 
gondola  axis  is  generally  not  aligned  with  the  true  vertical. 

The  vector  sum  of  the  three  accelerometers  is  in  a  direction,  and  has  a 
magnitude  that  must  be  the  vector  sum  of  the  gravitational  acceleration  (true 
vertical)  and  the  instantaneous  kinematic  acceleration.  Since  the  gondola  cannot 
turn  upside  down,  and  it  does  act  as  a  pendulum  which  has  a  restoring  force  in 
the  direction  of  the  gravity  vertical,  then  to  first  approximation,  the  average 
vector  sum  of  the  accelerometers  is  in  the  direction  of  the  true  vertical.  Applying 
a  forcing  function  (e.g.,  wind,  atmospheric  density  variations)  to  the  balloon 
system  will  alter  the  measured  accelerations,  but  unless  the  forcing  function  is 
applied  for  a  period  of  time  comparable  to  the  averaging  period,  the  average 
direction  of  kinetic  acceleration  must  still  be  zero.  The  horizontal  ucceierulions 
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are  immediately  countered  with  a  restoring  force  from  pendulous  swing,  while  the 
vertical  accelerations  are  countered  when  the  balloon  achieves  its  new  float 
altitude.  Thus,  we  expect  that  over  a  few  seconds  or  more,  the  horizontal 
accelerations  will  average  to  zero,  but  the  vertical  may  not.  In  both  cases,  the 
sensed  accelerations  still  average  to  the  direction  of  t'»p  true  vertical.  Affecting 
both  direction  and  magnitude  are  axis  misalignments,  biases,  and  drifts  of  the 
individual  sensors.  All  these  considerations  apply  for  translational  motion  only; 
they  are  significantly  more  complicated  when  rotation  aspects  are  included.  Any 
rotation  adds  centrifugal  or  centripetal  acceleration  in  a  variety  of  directions,  and 
coning  (i.e.,  the  motion  of  a  given  accelerometer  axis  describing  a  cone)  is  a  real 
concern. 

Combining  all  of  the  above  effects,  the  vector  sum  of  the  accelerometers 
generally  points  in  the  direction  of  the  gravity  vertical,  but  the  components  of 
that  vector  are  contaminated  by  axis  misalignments,  biases,  and  rotational  motions. 
Any  differences  detected  between  the  Z-accelerometer  vector  and  the  accelerometer 
vector  sum  must  be  attributed  to  a  combination  of  all  these  effects. 

The  one  identified  kinematic  acceleration  that  does  not  average  to  zero  is  due 
to  pendulum  motion.  Such  swinging  adds  an  average  acceleration  component  along 
the  pendulum  axis  (where  the  measurements  are  done  in  a  strapdown  system),  thus 
adding  to  the  apparent  measured  acceleration  in  the  direction  of  the  true  vertical. 
The  equations  describing  the  average  acceleration  are; 

ar(avg  =  g  [3  sin  0Q/0O  -  2  cos  0Q]  (4.26) 

a„(avg)  =  g  t  3  -  (cose  / 0  )(2  In  tan(P/4  +  0/2))]  (4.27) 

where  ar  and  az  are  the  radial  (along  the  pendulum  line)  and  vertical  accelerations, 
respectively,  g  is  the  ambient  gravitational  acceleration,  and  0Q  is  the  maximum 
swing  angle.  Figure  4.11  plots  the  ar  and  az  for  a  maximum  swing  angle  0Q. 

A  comparison  of  the  individual  accelerometer  signals,  then  transformed  into 
spherical  coordinates,  is  shown  in  Figure  4.12,  with  a  smoothed  version  shown  in 
Figure  4.13.  These  data  are  for  the  102.4  sec  period  during  the  quiet  time. 

During  this  time,  the  vector  sum  shows  a  total  acceleration  varying  over  about  500 
mGals.  The  Z-axis  misalignment  is  about  0.2  deg  with  a  variation  of  less  than 
0.01  deg.  Filtering  with  the  0.1  Hz  low-pass  filter  lowers  the  apparent  variation 
to  0.001  deg.  The  simplest  way  to  look  at  this  is  to  think  of  an  overall  Z-axis 
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Figure 


misalignment  of  0.2  deg  with  a  pendulous  swing  of  0.01  deg.  The  error  due  to 
the  misalignment  is: 

Ag  =  g  (1  -  cos  0^  (4.28) 

where  Ag  is  the  error  in  mGals,  g  is  the  gravitational  acceleration  in  mGals  and  0^ 
is  the  alignment-error  angle.  From  Figures "4.11  and  4.12,  and  Eq.  4.28,  the 
latter  introduces  an  error  of  much  less  than  1  mGal,  while  the  misalignment 
introduces  an  error  of  -6  mGals,  The  sum  of  the  average  errors  due  to  misalign¬ 
ment  and  pendulum  swinging  is  then  about  -6  mGals. 

Rate  Gyros 

The  gondola  gyros  also  form  an  orthogonal  system.  To  a  first  approximation, 
the  Z-gyro  shows  rotation  about  the  gondola  axis,  including  the  effect  of  the 
coupling  of  angular  momentum  between  the  balloon  and  gondola.  To  be  more 
accurate,  it  is  necessary  to  convert  the  gyro  data  into  an  inertial  coordinate 
system.  The  proper  approach  is  through  the  use  of  Euler  angles.  Although  some 
progress  has  been  made  in  this  area,  for  this  report  we  do  not  consider  Euler 
angles,  but  look  at  the  raw  data  alone.  From  a  practical  point  of  view,  the  X-axis 
and  the  Y-axis  gyros  must  average  to  zero,  as  the  gondola  cannot  turn  upside  down, 
or  even  depart  significantly  from  vertical  for  any  significant  amount  of  time.  This 
is  not  what  we  see  in  the  data  (Figures  4.14  and  4.15),  however.  What  the  X- 
and  Y-gyros  do  show  is  a  bias  which  must  be  incorrect;  therefore,  the  data  can 
reasonably  be  detrended.  Given  that  two  gyros  have  biases,  we  must  assume  that 
the  Z-gyro  is  also  biased.  The  magnetometers  and  the  compass  can  be  used  to 
calibrate  the  Z-gyi'o  by  computing  the  angle  to  magnetic  North  and  differentiating. 
Thus,  in  principle,  we  have  a  way  to  compensate  for  the  biases  in  all  three  gyros. 

Magnetometers 

Finally,  the  gondola  package  contains  a  three-axis  magnetometer  and  a  com¬ 
pass.  These  instruments  can  be  used  to  gain  information  on  rotational  aspects  of 
the  gondola  since,  given  a  good  earth's  magnetic  field  model,  the  data  will  yield 
absolute  angle.  Mathematical  continuation  models  also  apply  to  magnetic  fields,  and 
the  magnetic  data  available  here  at  30  km  altitude  is  unique  and  may  form  a  rather 
unique  data  set.  These  instruments  give  a  complementary  data  set  to  the  rate 
gyros,  and  allow  for  consistency  checks  and  noise  analysis. 

Two  sets  of  data  are  presented  to  show  three-axis  rotation  of  the  gondola: 
the  compass  data  showing  rotation  about  the  gondola  Z-axis,  and  the  magnetometers 
showing  all  three  axes.  Figure  4.16  shows  the  angle  between  the  gondola  X-axis 
and  magnetic  North  for  the  compass,  and  between  the  X-axis  and  true  North  for 
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figure  4.15  Example  of  Smoothed  3-Axis  Rate  Gyro  Data  and  Detrended  Version 


Figure  4.16  Comparison  of  Raw  Compass  Data  (to  Magnetic  North)  With  Raw  Magnetometer 
Estimate  of  Compass  Angle  (to  True  North).  Magnetic  North  is  about  +10  deg  to  true  North 


the  magnetometers.  The  magnetometer  should  show  about  10-deg  less  than  the 
compass  due  to  the  magnetic  azimuth.  There  is  still  a  30-deg  difference  between 
the  two  instruments.  The  reason  for  this  offset  has  not  been  determined  yet.  We 
have  decided  not  to  work  this  problem  for  this  flight,  at  this  time,  but  may  return 
to  it  depending  on  .he  results  of  DUCKY  II,  where  we  expect  much  higher  data 
quality.  Figure  4.17  shows  the  residual  magnetic  data  when  the  NASA  MGST-481 
model  is  subtracted  (as  based  on  the  WSMR  radar  position  data).  The  radial 
component  shows  a  400  nT  variation.  The  NASA  field  model  is  only  good  to  about 
1000  km  resolution,  and  does  not  account  for  local  effects  (e.g. ,  ionospheric  cur¬ 
rents,  solar  activity).  In  addition,  the  magnetometers  may  not  be  calibrated  to 
this  accuracy.  The  9  component  will  be  described  in  more  detail  in  the  next 
section.  The  $  component  shows  a  vertical  misalignment  of  about  1.35  deg  with  a 
variation  of  about  0.01  deg.  The  misalignment  does  not  agree  with  the  accelerom¬ 
eters,  but  the  short-term  variations  do.  Again,  the  magnetometer  data  for  this 
component  is  contaminated  by  the  same  effects  as  for  the  radial  component.  Clearly, 
the  magnetometer  shows  much  promise,  scientifically  on  its  own  accord,  as  well  as 
for  gondola-motion  determination.  Much  remains  to  be  resolved,  however,  and  for 
now  the  problems  will  be  shelved.  In  any  case,  the  data  set  is  unique  scientifi¬ 
cally,  and  should  be  looked  at  again. 

Z-Rotation 

Three  independent  sets  of  data  are  available  to  show  the  rotation  of  the 
gondola  about  the  Z-axis:  1)  compass,  2)  magnetometer  and  3)  Z-gyro.  Figure 
4.18  shows  the  three  sensors,  with  the  integral  of  the  Z-gyro  instead  of  the  rate 
data.  The  integrated  Z-gyro  data  has  no  calibration  point  as  this  must  come  from 
either  the  compass  or  magnetometer,  which  do  not  agree.  So  the  first  integrated 
Z-gyro  data  point  was  defined  as  zero.  Subtracting  the  three  data  sets  in  the 
three  combinations,  with  an  arbitrary  zero  point  is  shown  in  Figure  4.19.  The 
best  agreement  comes  from  the  magnetometer  and  integrated  Z-gyro,  at  least  for 
the  short-term  variations.  We  must  conclude  that  the  calibration  of  the  compass 
data  is  suspect,  but  all  differences  and  errors  in  rotation  data  are  suspect  as 
well.  These  data  clearly  need  close  scrutiny,  but,  again,  have  been  shelved  for 
now  as  the  sensors  for  DUCKY  II  have  been  improved  significantly. 

4.5.3  SUMMARY 

In  practice,  the  noise  evident  in  all  instruments,  and  the  need  to  convert  the 
data  into  an  inertial  frame,  make  this  a  difficult  problem.  Due  to  the  difficulties 
and  dead  ends  encountered  in  the  previously  described  portions  of  this  work 
(i.e.,  building  a  data  set,  analysis  of  translational  motions),  the  rotational  aspects 
have  not  been  rigorously  analyzed.  This  section  does  make  the  attempt  to  show 
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what  insights  have  been  gained  so  far  in  the  rotational  motions  of  the  balloon- 
gondola  system;  and  especially  to  summarize  what  we  perceive  to  be  necessary  in 
subsequent  flights.  Sufficient  information  was  derived  from  this  qualitative  and 
semi-qualitative  analysis  to  affect  the  design  of  DUCKY  II.  Suggestions  include: 

1)  Adding  a  swivel  (rotating  about  the  Z-axis)  between  the  parachute  and  gondola. 
This  should  help  uncouple  rotation  of  the  balloon  (presumably  driven  by  winds), 
and  gondola  rotation.  2)  Leaving  the  strapdown  system  as  is  for  the  foreseeable 
future.  The  gondola  is  stable  enough  to  consider  it  a  stabilized  platform  to  0.01 
deg.  3)  Improving  and  calibrating  the  compass,  magnetometers,  and  gyros. 


5.  CONCLUSIONS 
5.1  Summary 

The  flight  of  DUCKY  la  was  a  complete  success.  All  sensors  worked  as 
expected,  and  the  data  analysis  has  progressed  well  and  has  given  valuable  infor¬ 
mation  for  the  design  of  future  flights.  This  project  was,  and  still  is,  a  feasi¬ 
bility  project  where  the  goal  is  to  design  an  operational  instrument  to  measure 
gravity  at  high  altitudes;  that  is,  at  about  30  km  (100,000  ft).  The  original  plan 
calls  for  two  to  four  flights  total  to  understand  the  experiment,  and  actually  carry 
out  gravity  measurements  to  an  accuracy  of  1  mGal.  For  this  flight,  with  a  data 
analysis  scheme  not  yet  fully  developed,  and  inadequate  tracking,  we  feel  we 
disagree  with  the  DMA  gravity  model  by  26.5  ±  2.5  mGals.  For  a  first  attempt, 
this  is  a  very  good  result. 

The  strategy  in  data  analysis  was  to  focus  attention  on  "quiet  times"  where 
balloon  and  gondola  motions  were  obviously  the  lowest.  If  necessary,  the  analysis 
could  encompass  the  entire  flight.  Just  the  analysis  of  the  "quiet  times,"  how¬ 
ever,  proved  sufficient  to  justify  continuing  the  experiment  and  help  redesign  the 
system  for  the  second  flight.  In  particular,  three  issues  proved  more  difficult 
than  anticipated:  1)  tracking,  2)  rotation  sensors,  and  3)  a  total  inertial  navi¬ 
gation  solution.  Many  approaches  were  tried  to  address  these  issues,  and  the 
following  conclusions  were  reached.  The  radar  tracking  data  was  inadequate,  and 
was  totally  dependent  on  flying  the  experiment  near  high-quality  radar  instal¬ 
lations.  It  seems  that  the  best  solution  is  to  convert  to  GPS  tracking,  and  this 
will  be  done  for  DUCKY  II.  This  will  allow  for  performing  the  experiment  any¬ 
where  on  earth,  independent  of  ground-based  tracking  support.  Preliminary 
indications  from  other  GPS  experimenters  suggest  an  order  of  magnitude  improve¬ 
ment  in  tracking  data  quality.  The  rotation  sensors  were  particularly  sensitive  to 
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bias  and  drift.  Since  we  expect  noise  due  to  rotation  to  be  in  the  tens  of  milli- 
gals,  as  this  was  our  final  error  for  this  flight,  closer  attention  must  be  paid  to 
this  issue.  For  the  next  flight,  the  sensitivity  and  calibrations  for  these  in¬ 
struments  will  be  impi-'ved  by  about  a  factor  of  four.  The  total  solution  required 
an  all-encompassing  inertial  navigation  solution,  where  tracking  and  motion  data 
are  known,  but  the  gravitational  acceleration  is  not.  This  is  a  variation  of  typical 
inertial  navigation  systems,  but  is  capable  of  being  done.  This  part  of  the  analy¬ 
sis  will  be  submitted  to  the  private  sector  for  competitive  bids. 

Correcting  for  translational  motion  only,  we  find  the  average  measured  gravi¬ 
tational  acceleration  is  about  30  to  35  mGals  lower  than  expected  from  the  models. 
Preliminary  analysis  of  rotation  data  suggests  six  of  those  milligals  can  be  ac¬ 
counted  for  by  VSA  misalignment  within  the  gondola,  and  the  gondola  axis  being 
misaligned  with  the  true  vertical.  The  misalignment  is  about  0.2  deg.  In  sum,  we 
find  our  measurements  are  lower  than  expected  from  gravity  models  by  about  24 
to  29  mGals.  We  expect  most  of  this  residual  error  lies  in  noisy  tracking  data; 
instrument  biases,  drifts  and  misalignments;  and  incomplete  data-analysis  software. 

An  interesting  concept  has  developed  during  this  analysis.  The  data  from 
these  flights  happen  to  fall  in  a  potentially  interesting  region  for  General  Rela¬ 
tivity  experiments.  It  is  possible  that  the  gravitational  constant  G  varies  with 
time  and/or  distance.  Many  measurements  and  tests  have  already  been  done,  with 
insignificant  results.  No  measurements  have  been  done  at  the  distances  involved, 
from  earth  to  balloon.  Being  the  first  to  measure  gravity  in  this  region,  we  have 
the  potential  to  contribute  to  General  Relativity  experiments  involving  G(r).  It  is 
clear,  however,  that  we  must  break  the  1  mGal  level,  and  be  able  to  do  this 
during  ascent  so  we  have  many  samples  at  varying  distance.  How  well  we  can  do 
this,  or  expect  to  do  this,  remains  to  be  seen. 

In  summary,  the  flight  and  data  analysis  look  good.  We  have  specific  improve¬ 
ments  for  subsequent  flights,  and  we  are  going  ahead  with  the  project. 
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LIST  OF  ACRONYMS 


AFGL 

ARL 

CIGTF 

CiRiS 

DMA 

FM 

GPS 

I*S*P 

MSL 

NSWC 

OSU 

PCM 

PSL 

UTIAS 

VSA 

WSMR 

WHOI 

Z-100 


Air  Force  Geophysics  Laboratory 

Applied  Research  Laboratories,  University  of  Texas 

Central  Inertial  Guidance  and  Testing  Facility 

Completely  Integrated  Range  Instrumentation  System 

Defense  Mapping  Agency 

Frequency  Modulation 

Global  Positioning  System 

Interactive  Signal  Processing  (Bedford  Research  Associates) 

Mean  Sea  Level 

Naval  Surface  Weapons  Center 

Ohio  State  University 

Pulse-Code  Modulation  \ 

Physical  Sciences  Laboratory,  New  Mexico  State  University 
University  of  Toronto  Institute  for  Aerospace  Studies 
Vibrating  String  Accelerometer  (gravimeter) 

White  Sands  Missile  Range 
Woods  Hole  Oceanographic  Institute 
Zenith  Z-100  microcomputer  system 


Appendix  A 

Relative-Gravity  Measurements  et  Holloman  AFB,  New  Mexico 

Relative-gravity  measurements  were  made  at  selected  sites  at  Holloman  AFB, 
NM11,  in  support  of  the  AFGL  balloon-borne  gravity  experiment.  The  site  locations, 
chosen  in  accordance  with  the  experiment  requirements,  were: 

1)  The  Advanced  Inertial  Test  Laboratory  (AITL),  Building  1256:  This  is 

12 

an  absolute  site  measured  by  AFGL/LWG  in  1980,  the  Instituto  di  Metrologia  "G. 

Colonnetti"  (IMGC)  Italy  in  1980,^  and  the  Joint  Institute  for  Laboratory  Astro- 

14 

physics  (JILA)  in  1982.  The  absolute  gravity  marker  in  Room  10  in  the  AITL 
building  was  the  starting  and  closing  site  for  all  measurements. 

2)  The  High  Bay,  Building  850:  The  experiment  package  assembly,  testing, 
and  instrument  calibration  were  performed  in  the  High  Bay. 

3)  The  Altitude  Test  Chamber,  Building  1261:  Calibration  data  were  taken 
here  during  high-altitude  simulation  tests. 

11.  Iliff,  Robert  L.,  and  Sands,  Roger  W.  (1984)  Relative  Gravity  Measurements 

at  Holloman  AFB,  New  Mexico,  AFGL  Technical  Memorandum  Series  95 . 

12.  Iliff,  R.,  and  Sands,  R.  (1983)  The  AFGL  Absolute  Gravity  Measuring 

System,  A  Final  Report  and  Operating/Maintenance  Manual , 

AFGL-TR-83-0297,  ADA147853. 

13.  Marson,  I.,  and  Alasia,  F.  (1980)  Absolute  Gravity  Measurements  in  the 

United  States  of  America,  AFGL-TR-81-O052,  ADA699017. 

14.  Zumberge,  M.  A.,  Fowler,  J.  E.,  and  Gschwind,  J.  (1983)  Results  from  an 
absolute  gravity  survey  in  the  United  States,  J.  Geophys.  Res.  88  (No.  B9). 


4)  Secondary  Test  Chamber,  Building  1261:  Due  to  a  malfunction  in  the 
main  test  chamber,  a  secondary  high-altitude  simulation  chamber  was  used  for  a 
portion  of  the  tests.  The  floor  of  the  secondary  chamber  is  about  50  cm  higher 
than  the  floor  of  the  main  chamber,  so  the  gravity  value  is  lower.  The  gravity 
value  in  the  secondary  chamber  was  calculated  using  the  previously  determined 
gradient  of  0.285  mGal/meter. 

5)  Nenninger  Site,  Launch  Pad:  Measurements  were  made  30  m  diagonally  in 
from  each  of  the  four  corners,  and  at  the  center  of  the  launch  pad.  Five  points 
were  made  primarily  because  of  the  relatively  large  gravity  differences,  approxi¬ 
mately  one  mGal,  encountered  at  the  launch  pad  and  the  uncertainty  of  the  launch 
point.  Due  to  wind  shifts,  the  launch  point  frequently  is  not  known  until  just 
before  launch. 

It  should  be  pointed  out  that  the  measurement  made  at  the  center  is  not  at 
the  crosspoint  of  the  E-W/N-S  reflecting  markers  embedded  in  the  launch  pad. 
These  markers  no  longer  cross  at  the  center,  due  to  additional  lengthening  of  the 
pad  in  the  north-south  direction. 

The  gravity  transfers  were  made  from  the-AITL  building  with  a  LaCoste- 
Romberg  G-120  meter  loaned  to  AFGL  by  the  Geodetic  Survey  Squadron  (GSS). 

The  G  meter  was  set  on  its  leveling  feet  since  no  disk  was  available.  Further¬ 
more,  the  meter  was  placed  in  an  east- west  orientation  with  the  two  leveling  screws 
to  the  east,  and  the  single  screw  to  the  west  side.  Three  readings  were  made  at 
each  site  and  the  results  were  corrected  for  earth  tides,  meter  factor,  and  meter 
drift.  Meter  drift  is  based  on  constant  drift  from  start  to  closure  of  readings. 

Relative  Measurements 

The  confidence  level  of  the  relative  measurements  is  +  .015mGal  (except  where 
noted). 


1) 

AITL  (Building  1256)  (absolute  gravity  site): 

979139.600 

+  .009  mGal 

2) 

High  Bay  (Building  850): 

979136.002 

mGal 

3) 

Test  Chamber  (Building  1261): 

979136.536 

mGal 

4) 

Secondary  Test  Chamber: 

979136.393 

+  0.050  mGal 

5) 

Nenninger  Site  (Launch  Pad): 

a.  Southeast  Corner: 

979140.708 

mGal 

b.  Northwest  Corner: 

979141.782 

mGal 

c.  Southwest  Corner: 

979141.132 

mGal 

d.  Northeast  Corner: 

979141.286 

mGal 

e.  Center  of  Pad: 

979141.236 

mGal 
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Appendix  B 

Gravity  Model  Verification 
PCM/FM  Telemetry  Data 


Link  No.  1 


2258.5  MHz 


Encoder: 

Bits  per  Word: 

Bit  Rate: 

Output  Code: 

Frame  Sync  Pattern: 
Words  per  Frame: 
Digital  Words: 


Conic  410-2  S/N  1.062 
10 

128  Kbits 
BIO 

11111010111100110010 
47  (data) 

3 


Word  No. 
001 
002 
003 
004 
005 
006 
007 
008 
009 


Function 


Digital  Word  1  (see  digital  bit  assignments) 
Digital  Word  2  (see  digital  bit  assignments) 
Digital  Word  3  (see  digital  bit  assignments) 
Not  Available 
Encoder  Temperature 
Primary  Control  Unit  Temperature 
Primary  Command  Receiver  Signal  Strength 
0  -  15  PSIA  Altitude  Xducer 
0-2  PSIA  Altitude  Xducer 
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Word  No. 


Function 


010 

Oil 

012 

013 

014 

015 

016 

017 

018 

019 

020 

021 

022 

023 

024 

025 

026 

027 

028 

029 

030 

031 

032 

033 

034 

035 

036 

037 


0-0.2  PSIA  Altitude  Xducer 
Primary  30-Volt  Power  Monitor  (X8) 

Primary  12- Volt  Power  Monitor  (X3) 
Accelerometer  Z  Axis  ±  7.5g 
Accelerometer  Z  Axis  ±  2.5g 
Accelerometer  Y  Axis  ±  2.5g 
Accelerometer  X  Axis  ±  2.5g 
O  Volts  Calibration 

Backup  Command  Receiver  Signal  Strength 
Backup  30-Volt  Power  Monitor 
Backup  System  Internal  Temperature 
Link  No.  1  XMTR  Temp. 

Link  No.  2  XMTR  Temp. 

AFGL  Battery  Section  Temp. 

UTIAS  Battery  Section  Temp. 

VS  A  Compartment  Temp. 

Digital  Magnetic  Compass  Heading 
Command  Verification  Backup  System 
Command  Verification  Primary  System 
Blank  Tilt  Sensor  #113  Port  and  Starboard 
Blank  Tilt  Sensor  #114  Forward  and  Aft 
UTIAS  X  Gyro  Low  Gain 
UTIAS  Y  Gyro  Low  Gain 
UTIAS  Z  Gyro  Low  Gain 
UTIAS  X  Accelerometer  Low  Gain 
UTIAS  Y  Accelerometer  Low  Gain 
UTIAS  Z  Accelerometer  Low  Gain 
UTIAS  X  Accelerometer  High  Gain 
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Function 


Wor.d  No. 


038 

UTIAS  Y  Accelerometer  High  Gain 

039 

UTIAS  ?.  Accelerometer  High  Gain 

040 

UTIAS  X  Magnetometer 

041 

UTIAS  Y  Magnetometer 

A40 

UTIAS  Z  Magnetometer 

043 

UTIAS  Oven  Temperature 

044 

UTIAS  Box  Temperature 

045 

UTIAS  X  Gyro  High  Gain 

046 

UTIAS  Y  Gyro  High  Gain 

047 

UTIAS  Z  Gyro  High  Gain 

048  through  062 

Not  Used 

Word  No. 


001 


DIGITAL  WORD  ASSIGNMENT 


Bit  No. 

Indicated  Function 

01 

Flight  Time  Remaining  Tens  -2 

02 

Flight  Time  Remaining  Tens  -2 

03 

Flight  Time  Remaining  Tens  -2 

04 

Flight  Time  Remaining  Tens  -2 

05 

Flight  Time  Remaining  Units  -2 

06 

Flight  Time  Remaining  Units  -2 

07 

Flight  Time  Remaining  Units  -2 

08 

Flight  Time  Remaining  Units  -2 

09 

Termination  by  Primary 

10 

Termination  by  Backup 

002 


003 


01  CMD19  POWER  CONTROL  INDICATOR 

02  CMD19  POWER  CONTROL  INDICATOR 

03  CMD15  POWER  CONTROL  INDICATOR 

04 
05 
06 
07 
08 
09 
10 


01  Ballast  Pour 

02  Not  Used 

03  Not  Used 

04  Helium  Valve  Open 

05  Burst  Switch  Armed 

06  Burst  Pin  Pulled  Above  lOKft 

07  Burst  Pin  Pulled  Above  lOKft 

08  Not  Used 

09  Not  Used 

10  Not  Used 


WHOI  PRIMARY 
WHOI  BACKUP 
UTIAS  PRIMARY 
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GRAVITY  MODEL  VERIFICATION 
PCM/FM  TELEMETRY  DATA 


Link  No.  2 

2279.5  MHz 

Encoder: 

Bits  per  Word: 

Bit  Rate: 

Output  Code: 

Frame  Sync  Pattern: 

Words  per  Frame: 

20K-32  (Lab  Built) 

10 

19.2  kHz  Fixed 

B10 

11101101111000100000 

30  -  27  Analog  -  3  Digital 

Word  No. 

Function 

001 

Digital  Word 

1 

002 

Digital  Word 

2 

003 

Digital  Word 

3 

004 

Plate  Temp 

005 

Heat  Sink  Temp 

006 

Outer  Can  Temp 

007 

Inner  Can  Temp 

008 

Outer  Can  Control  Voltage 

009 

Inner  Can  Control  Voltage 

010 

Outer  Can  Power 

on 

Inner  Can  Power 

012 

Battery  Voltage 

013 

Chamber  Temp 

All  following  words  are  repeats  of  above  words 


on - 

Cardcage  Temp 

015 

Outer  Can  Temp 

016 

Inner  Can  Temp 

017 

Outer  Can  Control  Voltage 

018 

Inner  Can  Control  Voltage 

019 

Outer  Can  Power 
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GRAVITY  MEASUREMENTS  FLIGHT  H83-19 
STRIP  CHART  DATA  RECORDING  ASSIGNMENTS 


Recorder  #1 

Link  #1 

Word 

Primary  RCVR  Sig  Strength 

007 

Backup  RCVR  Sig  Strength 

018 

Command  Verify  Primary 

028 

Command  Verify  Backup 

027 

Digital  Mag  Compass 

026 

Z  Accel  ±  7.5g 

013 

Z  Accel  ±  2.5g 

Time  Code 

014 

Recorder  #2 

Link  #1 

Word 

Z  Accel  Low  Gain 

036 

Z  Accel  High  Gain 

039 

X  Accel  Low  Gain 

034 

X  Accel  High  Gain 

037 

Y  Accel  Low  Gain 

035 

Y  Accel  High  Gain 

Time  Code 

038 

kecorder  #3 

Link  #2 

Word 

Flat  Plate  Temp 

004 

Heat  Sink  Temp 

005 

Outer  Can  Temp 

006 

Outer  Can  Control  Volts 

008 

Outer  Can  Power 

010 

Inner  Can  Temp 

007 

Inner  Can  Control  Volts 

009 

Inner  Can  Power 

Oil 

Recorder  #4 

Link  #1 

Word 

X  Gyro  Low  Gain 

031 

X  Gyro  High  Gain 

045 

Y  Gyro  Low  Gain 

032 

Y  Gyro  High  Gain 

046 

Z  Gyro  Low  Gain 

033 

Z  Gyro  High  Gain 

Time  Code 

047 

Recorder  #5 

Link  #1 

Word 

Digital  Mag  Compass 

026 

Tilt  Port/Starboard 

029 

Tilt  Forward/Aft 

030 

Z  Magnetometer  UTIAS 

042 

Z  Accel  +  2.5g 

Time  Code 

014 
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Appendix  C 

Gravity  Preload  Preflight  Schedule 


IN  HIGH  BAY  (NIGHT  OF  LAUNCH) 

TIME  completed 

UTIAS  install  ice  pack 

30  min 

0030 

WHOI  install  front  plate  VSA  section 

10 

0040 

AFGL  install  batteries 

20 

0100 

(Note:  parachute  will  be  in  flight  trim  and  final  rigging  at  the  triplate) 

AFGL  run  preliminary  flight  line  checks 
(Refer  to  High  Bay  Checklist) 

30 

0130 

AFGL  install  1-in.  insulators,  web 
belts,  etc. 

10 

0140 

WHOI  perform  level  tests 

30 

0210 

AFGL  complete  High  Bay  checks 

AFGL  ensure  VSA  and  battery  sections  doors  in  place 
AFGL  power  down  except  for  WHOI 
AFGL  remove  auxiliary  battery  connector 
AFGL  secure  turnbuckles  with  safety 

wire  30  0240 

AFGL  install  all  doors  and  insulation  with  straps 


R 


TURN  OVER  SYSTEM  TO  CREW  CHIEF  FOR  WEIGHOFF  AND  TRANSPORT  TO 
NENNINGER  SITE 

60  0340 

AT  LAUNCH  SITE 

AFGL  remove  battery  compartment  door 
WSMR  perform  transponder  tests 

AFGL  do  control  center  tests  (Refer  to  Flight  Line  Checklist) 

AFGL  get  data  report  from  TM  re  sensors  status 

WHOI  record  data  for  15  min  60  0440 

Prepare  for  final  rigging  of  parachute  to  balloon 

AFGL  turn  power  off  manually  using  package  start  box 

Launch  crew  lay  out  balloon  and  parachute 

Connect  parachute  to  balloon  and  check  payload  release  device  and  burst  switch 
installation 

Ch-'ck  all  connections  among  gondola,  parachute,  and  balloon 

Reapply  power  to  system 

Remove  WHOI  battery  pack 

Remove  start  box 

Close  door 


Insulate  gondola 

45  min 

0525 

Final  weighoff  scale  check 

15 

0540 

PERMISSION  TO  LAUNCH 

0600 
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$ 


Time  (GMT) 

0515 

0700 

0901 

0944 

0955 

1003 

1013 

1035 

1100 

1309 


1320 

1341 

1400 

1409 

1411.29 

1417 

1418 
1431 
1452 
1504 
1528 
1555 
1604 
1609 


Appendix  D 

Control  Center  Log 


H83-19 

Gravity 


WX  is  go.  Mission  is  go.  Recording  updated. 

Control  Center  opened. 

Launch  Crew  departing  for  Nenninger  Site. 

Launch  Crew  at  Nenninger  Site. 

Laying  out  Bubble. 

Inflation  pressure  confirmed. 

Beginning  balloon  layout. 

Balloon  layout  complete. 

Holding  for  winds. 

Trouble  with  radios;  all  parties  requested  to  go  to  F2  and  stay 
there. 

Still  holding  -  winds. 

Starting  inflation. 

Approximately  15  minutes  to  launch. 

Clearance  to  launch. 

LAUNCH. 

CR/R  and  three  radars  OJ  this  time. 

Balloon  at  lOKft. 

Balloon  at  20Kft. 

Balloon  at  40Kft;  heading  SE  at  112  knots. 

FH13  called  in;  will  call  back  at  1615. 

Ea’loon  at  70Kft. 

FH60  heading  to  Hobbs. 

Balloon  at  float. 

F13  heading  at  Hobbs  Airport.  HF  check  showed  5707.0kHz 
worked.  Heard  60  weak,  but  could  understand  what  they  said. 
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1708 

1804 

1805 
1826 
1828 
1857 

1907 

1908 
1921 
1929 
1933 
2000 
2009 
2017 
2025 
2030 
2032 
2037 
2040 

2045 

2046 

2051 

2054 

2059 

2106 

2107 

2116 

2119 

2124 

2146 

2311 


FH13  called,  standing  by  at  Airport;  PH  393-9275 
BS  disarmed. 

Valve  opened. 

FH13  called,  standing  by.  FH60  will  call  in  i  hour,  1900Z. 

Valve  closed. 

FH13  called;  FH60  will  get  airborne  and  stay  in  Hobbs  vicinity; 
FH  13  will  stand  by. 

Balloon  at  85.8KFT  and  dropping  slowly. 

Pouring  90  sec  of  ballast. 

Pouring  45  see  of  ballast. 

Pouring  20  sec  of  ballast. 

Termination  NOTAM  given  (-1  hour). 

FH60  in  the  air. 

Pouring  20  sec  of  ballast. 

FH13  heading  to  west  side  of  Hobbs  on  Hwy  180/62. 

5  min  to  termination. 

TERMINATION! 

Balloon  past  chute. 

Payload  at  6 OK. 

Payload  at  55K. 

Payload  at  45K. 

Balloon  impact;  080°  51.5  nm  Earlstrad  VOR. 

100°  57  nm  Roswell  VOR. 

Payload  at  35K,  loss  of  radar. 

FH60  lost  payload  in  clouds. 

Payload  at  25K, 

Pouring  ballast. 

TM  loss. 

Pour  command  dropped. 

13  and  60  have  a  visual,  payload  at  approx  7K. 

Payload  Impact:  084°  11.5  nm  Hobbs  VOR. 

300°  57  nm  Midland  VOR. 

Recovery  at  site. 

Picked  up  payload,  heading  to  Airport  to  drop.  Lt.  DePiero  off 
with  60,  then  picked  up  balloon. 

CLOSING  DOWN. 
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